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subroutine and high resolution precision plotting software 

to provide the user with a highly accurate time response of 

a system which has been designed on the OPTSYSX Frogran. 

Lis series of programs permits the user to rapidly design, 

analyze and test all types of Optimal Systems Control frob- Br 
lems. Examples of the various types of problems are worked oe 
through to illustrate all of the capabilities available. 
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state (Ns,Ns) or output (No,No) weighting matrix 
control (Nc,Nc) weighting matrix 

control gain matrix (Nc,Ns) 

centrol (No,Nc) or noise (No,Ng) feedforward 
matrix 

open-loop dynamics matrix (Ns,Ns) 

control distribution matrix (Ns, Nc) 

state disturbance distribution matrix (Ns,Ng) 
measurement scaling matrix (No,Ns) 

estimator gain matrix (Ns,No) 

number of ccntrols 

number of process noise sources 

number of states 

humber of okservations or measurements ie 
white process noise covariance matrix (Ng, Ng) 

white meas. noise covariance matrix (No,No) 

steady-state covariance matrix of control (Nc,Nc) 

control vector (Nc,1) er 
control input (Nc,1) . 
state vector (Ns, 1) 

state vector derivative (Ns,1) 

estimate of state vector (Ns,1) 

derivative of estimate of state vector (Ns,1) 

state reconstruction error (Ns,1) 


output/measurement vector (No, 1) 
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I. INTRODUCTION 


The rpurpose of this thesis is to describe and denon- 
strate the modification and additions to the existing 
FOETRAN proyram (OPTSYSX) which is used in the study, 
design, and application of Optimal Systems Control theory. 

The CPTimal SYStems control program (GPTSYS) was origi- 
nated in 1971 by Hall [Ref. 1] to support his research in 
rotary-wing aircraft control systems. The most recent 
program modifications were made by Walker [Ref. 2] and Liu 
[Ref. 3] of Stanford University and are designated OPTSYS 4 
and CFISYS 5, respectively. The OPTSYS modifications made 
ty Hoden [Ref. 4] were primarily devoted to creating a 
user~friendly interactive version (OPTSYSX) of the CPTSYS 4 
Frojram. 

The goal of this thesis work was to develop a program 
set which will operate in an interactive mode and ;lot the 
time resronse of a State Variable Control System which has 
been develoyed and/cr analyzed using the OPTSYSX Progran. 
Minor modifications tc the OPTSYSX Program were aecessary to 
allow the user to build a data fiie of the matrices required 
for the time response calculations plus additional matrices 
which cculd be used again upon reentry to the oOPISYSX 
Program without the laborious and time-consuming task of 
reentering each data element in all of the required 
matrices. 

It is assumed that the reader/user is familiar with the 
Tasic cencerts of Control Theory and Optimal Systems Design. 
The symkcl/naming conventions of Bryson [Ref. 5] are used in 
the projram operation discussion and in the exanuples of the 
yroblems solved using this systen. A glossary of the 
symbols and abbreviations used in this discussion is 


provided on page 8. 
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An explanation of the OPTSYSX capabilities and a program 


set overview are presented first. 

This work concludes with examples of variois types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final results. A set 
of complete program listings are included in Appendices A, 
E, C and D. 


11 


wate ae 


Sonn Shen... 


we 


A. GENERAL BACKGROUND 


OFTSYSX is an interactive, double-precision FORTRAN 
program employing modern control theory analysis techniques. 
Its extensive capabilities include the synthesis and anal- 
ysis of filters and regulators as well as eigensystem anal- 
ysis, modal distribution, transfer function analysis and 
power spectral density computations. The modifications to 
the OPTSYS Program introduced by this thesis work have not 
affected any of the program's original capakilities. 


B. PROGRAM OVERVIEW 


CETSYSX is an extremely large and complex program with 


over 2800 lines of code. In order to use this program in 
its small (set up for 32X32 matrices) version, the user is 
reguired to extend the IBM 370/3033 virtual machine (VM) 
memory Capacity beyond 720 kilobytes whic«c is th= default VM 
memory size. A significant increase in the size of the 
OPTSYSX program would make the program too large to operate 
On a che megabyte VM, the largest virtual memory available 
or. a user's virtual machine. The high resolution plotting 
software is limited to Singie pre sion Variables. 
Therefore dcukle precision liitrazy routines cannot be called 
from the plotting progran. For these reasons the task of 
obtaining the time response of a system was divided between 
three programs, OPTSYSX, OPTCALC and OPTPLOT. An Executive 
program (CPTSYS EXEC) was written to make the interfacing of 
the three programs transparent to the user. 

Minor modifications were made to the OPTSYSK progran 
including the additicn of three subroutines to handle the 


iuput and output of matrix data to and from a data file on 
the user's disk. The OPTCALC program performs the double 
precision numerical integration of the system »f equations 
over time and creates another disk data file of the state 
variakle variation with time. OPTPLOT takes this time 
response data and presents it in a graphical format on the 
TEK 618 graphical display or as a VERSATEC pen pilot. 


C. OFTSYS EXEC 


The CPTSYS EXEC is written in the EXEC 2 lanjuage. This 
language allows the EXEC to issue almost any command that 
can ke entered in the direct mode at the terminal. Therefore 
ar EXEC is the ideal controller for the "black box" type of 
system where the user is not aware of what is actually 
taking place within the program(s). The OPTSYS EXEC was 
written to complete all of the reguired interfacing between 
the three programs (OPTSYSX, OPTCALC and OPTPLOT), without 
the direct guidance or control of the user. By answering 
questions presented on the terminal screen, the user deter- 
mines the logic flow through the EXEC while the EXEC estab- 
lishes the appropriate FILFDEFs and loads the frograms 
Teguired by the user's desires. 


De. OFTSYSX MODIFICATIONS 


Three subroutines (RDMATF, RDMAT, and WRTMAT) were added 
to the CPTSYSX Program for data file read error check, 
matrix input from a data file and matrix output to a data 
file, respectively. These three subroutines provide the 
user with the opportunity to save the [F], [G], {H). 
[GAMMA], , A} and [BJ] matrices for use in a subse,uent run of 
tne CETSYSX Progran. The WRIMAT subroutine also saves the 
{C] and [(K] matrices for system time response calculation 
and ;lotting by the OPTCALS and OPTPLOT programs. 
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The PDMATF subroutine is used to check for the exis- 


tence of a previously generated file containing matrix data. 
Seven flags may be set by this subroutine. Six of these 
flags correlate with the six matrices that the user may save 
for reuse later in the OPTSYSX Progran. The remaining flag 
(IRDMAT) must be set to enable the RDMAT subroutine to read 
matrix information from the data file. aA READ statement of 
the form 


READ (9,111,EKR=222, END=333) A,3 


(where 111", "222" and "333" designate line numsers for the 
FORMAT statement and branch on ERROR or branch on 
END-CF-FILE routines, respectively) is used for the data 
file check. The nenexistence of the rile or premature 
END-OF-FILE are detected by the ERR and END shecks which 
cause a Tranch to a routine that sets the IRDMAT flag to "0" 
and returns to the calling progran. 

If no error is detected during the initial read 
attempt, the variable B is checked for the sentinel "1", 
This second check is to help ensure that the file is actu- 
ally a file which contains valid matrix elements. The user 
is then presented a message which asks if he/she wants to 
use the matrices which are available. The user may respond 
with on cf three answers: 


1. Use all of the matrices. 
2. Use selected matrices. 
3. Use none of the matrices. 


If the answer is "1" or "2, the subroutine reads the matrix 
dimensions (Ns, Nc, No andNg) from the data file and 
changes the IRDMAT flag to 1 to key the RDMAT subroutine to 
read the matrix elements from the data file. If the answer 
is "1" all of the matrix-save flags are set to "1", 


Vy 


If the answer is "2", the user is given the opportu- 


nity to select individual matrices for reuse while rejecting 
other matrix information. This is accomplished by setting 
individual matrix-save flags to "1" if the matrix is to be 
saved and "0" if new matrix data will be inpat from the 
terminal. If the answer is "3" (Use none of the matrices) 
the IRDMAT flag is set to 0 and the subroutine returns to 
the main program. When all actions have been conpleted, the 
flag infcrmation and the matrix dimensions are passed to the 
main ysrogram for later use. 


2- RDMAT Subroutine 


ee 


The RDMAT subroutine is used to read all of the 
matrix information in the data file and transfer the infor- 
Mation to the aprfropriate variables. As previously 
discussed, The acticns of this subroutine depend on the 
status of the IRDMAT flag. If this flag had been’ set to 
"0", no read operations are attempted and program flow iame- 
diately returns to the calling progran. 

When the IRIMAT flag is set to "1", the RDMAT 
subroutine reads the matrix dimensions from the data file, 
and uses these dimensions to transfer the matrix information 
from the file to the appropriate variables. The file matrix 
dimensions are used for the read operations and are not fed 
back to the calling rrogram, since the dimensions of some of 
the matrices which are not being reused may have changed 
from the the previous run. Similarly, using the current 
matrix dimensions in the RDMAT subroutine would cause iata 
read-in protlems due to the changing number of elements in 
€ach matrix as the matrix dimensions vary. 


3. WHIMAT Subroutine 


The WRTMAT sukroutine is used to write a data file 


of the djata file flags, the matrix dimensions ind selected 
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BRatrices. When the user has completed the analysis/design of 
the system of interest, the WRTMAT subroutine asks the user 
if heyshe wants to calculate the time response of the system 
which the user just designed. If the user answers YES, the 
WETMAT subroutine generates a data file of appropriate 
matrix information and halts execution of the oOPTSYSX 
program. Control then reverts to the OPTSYS EXEC. If the 
user answers NO, the WRTMAT subroutine returns control to 
the main program and normal OPTSYSX program oferation 
continues. 

The informaticn written to the data file consists of 
2"1%s (which are used as a sentinel or flag by the RDMATF 
subroutine (as previously explained) and in a similar manner 
ty the OPTSYS EXEC), followed by the matrix dimensions (Ns, 
Ne, Ne and Ng) and then by the [FJ], [G], [H], [GAMMA], [C], 
[Xk], {A] and [B] matrices. These matrix elements are written 
to the OPIMAT DATA file using a 4D20.13 format as a compro- 
mise Letween the maximum feasibie accuracy of data exchange 
ketween the double-precision programs and the use of a 
moderate amount of the user's disk space. 


E. OFTCALC PROG RAS 


The OPTCALC program is a FORTRAN interactive 
double-precision system integration routine. This program 
uses the International Mathematical 6& Statistical Library 

IMSL) subroutine DGEAR to perform the numerical integration 
of the system under analysis. The stiff system mode of DGEAR 
is used in order tc provide the capability to do time 
response calculations of the X-29A longitudinal axis back-up 


mode system which is an 98 X 98 stiff systen. 


werd 


2. System £ 


guation Representation 


The OPTCALC program uses the state variable format such as 


xdot = [F]*x + [3 }¥uc (2.1) 


to define the system. In this system the [F] matrix is the 
open~loog dynamics matrix (system or plant) and the [G] 
matrix is the control matrix. The variaole assignments are 
X as the state vector and uc as the control input vector. It 
follows that xdot is the time derivative of x. 

Various forms of equation 2.1 are used for all the 
tine resronse calculations. The [F] matrix is modified to 
{F+3*C]1 for closed-loop (regulator only) system calcula- 
tions as in ejuations 2.2 and 2.3. 


xdot = [ F+G#c J*x + (| G]*uc 


u = (C ]}*¥x + uc (2.3) 


For this closed loop system, the (C] matrix is the ccntrol 
gain cr regulator gain matrix and uis the total input 
vector. 

The combined filter and regulator systems can be 
represented by equations 2.2, 2.4, 2.5, 2.6 and 2.7. The 
(H] matrix is the measurement scaling (observer output) 
matrix and the [K] matrix is the estimator or Kalman filter 


cbserver gain matrix. The variables xe, xedot and X are the 


state estimate vector, the derivative of the state estimate 


vector and the state reconstruction error, respectively. 


1The OPTSYSX sign convention for the C matrix is the 
negative of the ~“standar normally _used in  centrols. 
Therefore PE+gee | has the correct sign for OPISY5X natrices. 
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xdot = [F+G4#c ]¥*x + [G]*uc (2.2) 


z = [H}*x (2.4) 


xedot = [FJ*#xe + [G]*u + [K]*(z - [H]*xe) (2.5) 


Cc 
" 


[C }*xe (2.6) 


% =x - xe (2.7) 
Ecuations 2.2, 2.4, 2.5, and 2.6 can be combined into the 
augmented matrix form of equation 2.8. 


j xdot | i F G*C 


1 lx | IG] (2.8) 
lxedot| ~ | xan Fecac-xen |" txe]*| 


*1C 
3 


Equation 2.8 is an augmented equation in which the Ns dimen- 
sion has been doubled andthe state and state estimate 
vectors have been combined into one vector of 2*Ns length. 
The [G] matrix is also augmented by repeating the first Ns 
rows again beginning at row Ns+1, making the new [G] matrix 
dimensions (2*Ns,Nc). 

In a Similar manner, a filter only systen can repre- 
sented ty tke same equations with the [C] matrix set to 0. 
These equations are: 


xdot = (FJ*x + (G]*uc (2.1) 
z= [H)*x (2.4) 
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xedot = [F]*xe + [G]*u + °K]*(z - [H )*x2) (2.5) 
% =x - xe (2.7) 
The cemkEination of eguations 2.1, 2.4, 2.5 and 2.7 into an 
augmented [FJ] matrix is similar to equation 2.8 but with 


fewer terms and the upper right yuadrant eyual to zero. 
This filter-only augmented system equation is shown kelow. 


+f" le 


xe| 


xdot 
xedot 


F 0 
K*H F-K*H 


ae (2.9) 


G 


ae System Selection 


When the CPTCALC program is run, the |{Fj, [GI 
and |Cj matrices (and JH] and |K{ matrices, if available) 
are presented on the terminal as a check and as a reminder 
of the characteristics of the system that has been passed 
from OPTSYSX to OPTCALC. The user is requested to select the 
type cf system response to be calculated. 


1. CFEN LCOP TIME RESPONSE. 

2- CLOSED LOOP TIME RESPONSE. 

3. CLOSED LOOP FILTER ONLY TIME RESPONSE. 

4. CICSED LOOP FILTER + REGULATOR TIME RESPONSE. 


Selection of 2, 3 or 4 forms the appropriate system matrix 
equations 2.2, 2.9 or 2.8, respectively and doukles the 
length of the JG|{ matrix, if required. 


bk. Defininy Calculation Limits and Inputs 


After the user determines tne type of system 


under study, the OPTCALC program prompts for the integration 


start and stop times and the number of data points desired. 


The user has some ccntrol over the tradeoff batween curve 
fidelity and computer time used by varying the number of 
data points calculated. Computer time use is normally a 
factor only on very large systems. If less than 200 points 
are calculated, the OPTPLOT program uses a curve smoothing 
function which may cause minor inaccuracies in the plotted 
curve but avoids the sharp peaks and irregulac appearance 
generated by plotting straight lines between an insuificient 
number of data points. When 200 or more points are calcu- 
lated no smoothing is done. The points are connested by very 
short straight lines which has the appearance of a smooth 
curve. 

Step or ramp functions are available as control 
inputs. Only one type of function can be used for each 
contrel, but the magnitude and start and stop times can be 
set as desired. 

The state and state estimate initial conditions 
can Ee set to any value by the user. The control inputs and 
non-zero initial conditions can be used simultaneously, if 
desired. Before the time response calculations begin, the 
user is given the opfortunity to make changes in any area of 
the system integration initial conditions that have been 
previcusly selected. 


3. System Time Response 


Equation 2.1 is evaluated directly in the open loop 
system response calculations. The FCN subroutine was 
written to evaluate the system of equations for the DGEAF 
IMSL subroutine. Each time the FCN subroutine is called by 
IGEAF, it updates the control inputs (uc) and then evaluates 
each state derivative by summing all the terms across that 
row of the [F]*x and [G]*uc matrices. The same FCN sukrou- 
tine is used for all system integrations. As 2xplained in 
the ;revious section, the [F] matrix is replaced by the 
{F+3*C] watrix for clcsed loos system response problems. 
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a. Systems With Filters 


Augmented eguations 2.9 and 2.8 are ised for the 
time response evaluation of systems with filters only ani 
systems with filters plus regulators, respectively. The 
augmented matrix is developed as a dummy matrix and is then 
inserted as the [F] matrix witn the dimensions doubled 
(2*NS). The [G] matrix is also augmented by repeating the 
first Ns rowS again beginning at row Nst1, maxing the new 
[GS] matrix dimensions (2*Ns,Nc). 

The augmented system can be evaluatad by simply 
doubling the old system row and column dimension (Ns) and 
cailing the DGEAR integration subroutine. Using this 
metaod, the existing FCN subroutine reguires no changes to 


evaluate the augmented systen. 
E. OPTCALC Output 


The OPTCALC program uses FILEDEF 8 for the data 
file outfut as well as FILEDEF 5 to read and write tec the 
terminal. The output data file contains tke foilowing 
discrete information: the matrix dimensions Ns and Nc, the 
augmented matrix dimensions, the wuumber of data points 
calculated and a flag to indicate that an augmented matrix 
was calculated. The (C]) matrix is passed to permit the 
calculation of u the total control input to the system. The 
final fortion of the data file is individual data foints of 
time, external control input (uc) and each state (x) and 
State estimate (xe). This data file provides all of the 
data reguired by the OPTPLOT program to make a smooth graph- 
ical resyonse curve. 


F. OFTPLOT PROGRAS 


The CPIFLOT program is a FORTRAN interactive flotting 
program using the Display Integrated Software System and 
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Plotting Language (DISSPLA) ky Integrated Software Systems 


Corporation. This program provides the user a high resolu- 
tion grapfhical display of the system's time response ani if 
desired will provide a VERSATEC pen-plot of the same graph. 


1. General Operation 


Flotting data is received from the projram OPTCALC 
via a data file on FILEDEF 8. The types of data provided in 
the file are discussed in the previous section. The program 


presents the user with a series of guestions to determine: 


The number of curves to P19t 

Select the type of variable for each curve. 

Select the variable subscript for each curve. 
Select the number of headings and contents of each. 


eWwWhw 


The program then plots the selected variables and frovides 
the user a graphical display on the Tektronics 618 (TEK 618) 
display. 

The following Main Menu is then presented which 
Provides the user with the major decision points of the 


program. 
1. BEGIN NEW GRAPH CF OTHER CONTROLS ,STATES, OR ESTIMATES 
2. REFLOT PREVIOUSLY SAVED GRAPH DATA. 
3. EDIT THE CURRENT GRAPH. 
4. FLCI REVISED GRAPH ON THE TEKE18. 
5. uaa AND/OR MAKE METAFILE OF THE CURVES 
PREVIOUSLY SAVED. 

. The purpose of each of the selections is self- 


explanatory, however the methods of their use may not be. If 
number 3 is selected, the user is then presented the 
following Edit Menu of items to make additions, deleticns or 
corrections to the curves that are plotted on the TEK 618 


screen: 
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1. CHANGE VARIABLES Of ADD A CURVE ON THE CURRENT PLOT. 
2. DELETE CURVE FRCM CURRENT PLOT. 
3. EDIT CURVE aT EEAD 
4. EDIT PAGE HEADING ( ie 
‘a 5. CHANGE THE Y-AXIS SCALE. 
6. CHANGE THE TIME AXIS SCALE. 
7. CHANGE THE PLOT SIZE. 
8. CHANGE THE LETIFRING HEIGHT. 
9. CHANGE POSITION OF THE LEGEND. 
10. EDITING COMPLETE. 


This extensive list cf modification capabiliti2s provides 
the user with the tocl to make almost any imaginable altera- 
tion to an existing flot. 
os Since the OPTELOT program receives the time response 
data from the OPTCALC program, item 6 of the Edit Menu 
cannot be used to expand the time scale beyond the time span 
previously calculated. Therefore the time axis change 
e feature can be used cnly to select a subset of the original 
data. 


2. WVERSATEC Pen-plot 


a The VERSATEC pen-plots are provided through the 
a2 DISSPOP portion of DISSPLA. In order to use the DISSPOP 
feature, a device independent plot file called a metafile 
o must te generated. 

Laid To view the graphical time resporse plot on the TEK 
‘ 618 terminal, the TEK618 cption of DISSFLA must be called 
| within the OPTPLOT plotting program. The graphical image 


data is ther sent to the TEK 618 display screen. 

A metafile is created when the COMPRS option of 
DISSFLA is called by the OPTPLOT plotting proyram. When the 
Flotting, pregram is executed with the COMPRS option, the 
graphical image data is sent toa metafile on the user's 
disk. The TEK618 cfption and COMPRS o;tion are mutually 
exclusive (only one can be active at a time), therefore 
graphical data cannot go to both the terminal screen and the 
metafiie, concurrently. As a further complication, the 


TEK61€ and COMPRS oftions cannot be used in alternating 
pattern, first to originate and edit each graph and then to 
add this graph to a metafile possibly containing several 
cther graphs. Therefore if more than one pen-plot per 
terminal session is desired, some type of capability must be 
provided to save the information reguired to reproduce a 
given graph. 

Khen the user attempts to leave the current plot 
(ie. selecting items 1, 2 9r 5 of the Main Menu) the progran 
asks the user to save the current graph data for later use 
in generating a retafile. This feature provides the cafa- 
bility to save any desired graph data in order to later make 
a metafile and obtain a pen-plot. When the COMPRS oftion is 
used (by selecting item 5 of the Main Menu), any number of 
graphs tay be added to the metafile up to the limit of 
available user disk space (provided graph data has been 
previously saved). After exiting the OPTPLOT program, the 
CPTSYS CXEC asks the user if he/she wants a hard copy of the 
metafile that had been generated during the session. If the 
user answers YES, the OPTSYS EXEC cails the DISSPOP EXEC 
with the VRSTEC opticn. 

When the user exits the DISSPOP EXEC, the CFTSYS 
EXEC gives the user the options to: 


Whoo 
e a e 


The option to run OFISYSX again allows the user to use all 
er part of the matrices that had been saved in tae data file 
without manually reentering each element. The OPTCALC 
option could be exercised if the user wants to use the sane 
system matrices again, but change the control input of 


initial conditions cr change the type of system (open, 
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closed, filter only cr filter plus regulator) that was eval- 
uated on the previous run. 
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This chapter contains several basic examples of the four 
types of problems which may be solved using OPTSY SX, OPTCALC 
and OPIPLOT under control of the OPTSYS EXEC. Included with 


these examples are copies of each recorded terminal session. 


A. OPEN-LOOP SYSTEM TIME RESPONSE 


The following ofen-loop system example was taken from 
[Ref. 6, pp 5.3 - 5.7]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case letters or numbers 
below each "2", 


record on 
BEGIN RECORDING OF TERMINAL SESSION = 
KR; 7=0.01/0.92 19:58:26 
optsys 
THE OFTSYS EXEC CONTROLS A TRIO OF PROGRAMS: 
1. OPTSYSX FORTRAN SYSTEM ANALYSTI Ss) 
2. OPTCALC FORTRAN CALCULATE TIME RESPONS as 
3. OFTPLOT FORTRAN DISSPLA PLOTTING ROUTINE) 
EACH PROGRAM PASSES INFORMATION TO THE NEXT PROGRAM 
THROUGH A DATA FILE WRITTEN TO THE USERS DISK. IN THIS 
CASE, THESE FILES ARE “OPTHAT DATA" AND “OPTPLOT DATA". 
THE SIZE OF THESE FILES VARY WITH THE SYSTEM OKDER, AND 
CAN USE ABOUT 20% OF TEE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 
TYPE "E™" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 


PEt SSS TEST SEL SS SST ESSE TS TTS ST TS TSS TT TS TS TTS STS TTS SS TES EST SS SS 


YCU HAVE A DATA FILE NAMED SOPTMAT DATA" ON YOUR A 


DISK THAT KAS FREVIOUSLY GENERATED BY THE OPTSYS 
FROGRAM AND CCNTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 

IF YCU WOULD LIKE TC USE THESE SAME MATRICES FOR 
THIS RUN, THE OPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROFRIATE TIME, 


IF YOU TYPE (Y) ES. 
ANY OTFER INPUT WILL RESULT IN THAT FILE BEING ERASED! 
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PES PSE ER ESTE TESS STS ESET SS TE SESS SSE SEL SESE STE SSS SPSS ESE LSE ELST SE F 
Y 


ae Re Re eae ee ae a oe eae ake ah ae ake ee a ae a a ok ee a ee a RK RK KK kk ak kK Kk 


DO YOU WANT THE NUMERICAL OUTPUT ae OPTSYSX IO GO 
TO YOUR TERMINAL 44 °8E* 3} OR TO A D(ISK) FILE? 


ae a eK ee oe he oe eee he he eae ee eh oe ee ee KK RK KK KK a eK kK & 


s 


cree eet SSeS St SSS StS SESS STS SSS SES TLE SESE SSE SITS SSS SS SSE SEL STS SS 
OUTPUT WIIL COME TO YOUR TERMINAL SCREEN. 
PS SS SSESS SSS SSeS StS SSS SSS TE STS SL SSS ET ETS ST SS TST TSS SSS ES STS tS SS 
LOPES PERI 
CPISYSX IS A_COMPIE ERACTIVE OPTIMAL SY STEMS 
POLLORING TYPES. MS CONTROL EQUATIONS? 
XDOT = {F}*X + {G}¥*U + {GAM} * (W+WO) 
MEASUREMENT EQUATION-- 
Z = {B}*X + {D}*W + V 
REGULATOR PERFORMANCE INDEX-- 
1/2 * INTEGRAL (¥ *{A}*Y + U *{B} *U) DT 
STATE FEEDBACK GAIN DEFINITION-- 
= -{C}*X 
DO YOU WISH TO CONTINUE? TYPZ "YES" CR "NOM, 
--DATA ENTRY-- 
S SPECIFICALLY D 


z ESIG 
INTERACT EVERAL 
ABLE TO THSERS: 


Cy 
u 


NED TC READ 
ALTERNATE 


poker) 
WHO 


YSX 

TA 

VAIL 

THE "FV," Gt") AND "GAMMA" MATRICES 

READ FROM SEPARATE DATA FILES. 

THE Ror "Gt, AND "GAMMA" MATRICES MAY BE 
TLY DEF{NED’ WITHIN SUBROUTINE "SETUP", 
THEK CASE, itis USER ese Ch Opt ate A COPY 
EF EXA 

XAM 
ISH 


;OGRAM "LISTING AND 
LES CONTAINED IN S/R NSETUD", } 


TO CONTINUE? TYPE “YFS" OR SHON, 


OH UU KRU YH 


Oo mya et! 
hood 


DO YOU WISH TO INPUT THE "F", "G", AND “GAMMA" 
MATRICES FECM SUBROUTINE "SETUP" JAW THE 
METHOD DESCRIBED ON THE PREVIOUS SCREEN? 


TYPE "YES" OR "NO". 


27 


ee ie Se ae Se BR i ia a dias OE Ea SAO Dele Sete ete Tees a ie eee ita ieee a ee 


GENERAL OPTSYSX OPTIONS: 


OFTION 1 -- SYSTEM ANALYSIS WITHOUT 
OPEN-LOOP EIGENSYSTEM CALCULATIONS. 
CFTION 2 -- SYSTEM ANALYSIS WITH OPEN-LOOP 
EIGENSYSTEM CALCULATIONS. 
OFTION 3 -- OPEN-LOOP EIGENSYSTEM FOUND 
AND PROGRAM TERMINATES. 
{"F"-MATRIX ENTRY FOLLOWS IMMEDIATELY. } 
OEFTION & -- MODAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
OR STEADY-STATE ANALYSIS. 
SELECT AN OPTION: 1,2,3, OR 4. 


[CO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
TYPE "YES" CK "NO", 


CPEN-LCOP TRANSFER FUNCTION OPTIONS: 


CPTION 1 
CFTION 2 
CPTION:..3 
CPTION 4 


NO OPEN-LOCP TRANSFER FUNCTID NS COMPUTED. 
FCLES, RESIDUES, AND ZEROS COMPUTED. 
CNLY POLES AND ZEROS COMPUTED. 

CNLY POLES AND RESIDUES COMPUTED. 


SELECT AN OPTION: 1, 2, 3, OR 4. 


NCISE TRANS FER 


OPTION 1 
CETION 2 
CPTION 3 
CPTION 4 

SELECT 


COMPENSATOR 
CPTION 1 
CETION 2 
CPTION 3 
CPTION 4 


TR 


{NOTE: 


FUNCTION OEFTIONS: 

NC NOISE TRANSFER FUNCTIONS SOMPUTED. 
FCLES, RESIDUES, AND ZEROS COMPUTED. 
CNLY POLES AND ZEROS COMPUTED. 

CNLY POLES AND RESIDUES COMPUTED. 
OFTION: 1, 2, 3, OR 4. 


ANSFER FUNCTION OPTIONS: 

NO COMP. TRANSFER FUNCTIONS TOMPUTED. 
FCLES, RESIDUES, AND ZEROS COMPUTED. 
CNLY POLES AND ZEROS COMPUTED. 

CNLY POLES AND RESIDUGS COMPUTED. 


A COMPENSATOR TRANSFER FUNCTION CAN BE 
CCMPUTED ONLY IF BOTH A FEGULATOR 
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Oe owe 


AND FILTER ARE SYNTHESIZED 
AND/OK INPUT.} 


SELECT AN OFTION: 1, 2, 3, OR 4. 


WILL A FEED-FORWARD DISTRIBUTION MATRIX 
{"D" - MATRIX} BE INPUT ? 


TYPE "YES" OR "NO", 


THIS POLYNOMIAL ARE PRINTED CUT AND 
OW MANY EXTRA ROOTS EXIST AT Z 0. 
-IE} IS CONSIDERED ZERO. 


* 
ALUF OF THIS EERRARETER {IE} IS 6. 


= 1.0E-6. 


FFERENT MARKOV CRITERIA, 
ALUE. 


HE LEFAULT VALUE, TYPE "0" {ZERO} 


ae 


mw AMD BIA mHN 
is] Pa be oh 


POWER SPECTRAL DENSITY {PSD} OPTI 1 3 


OPTICN 1 -- TE ane pa OF TH 
LS OF TH 


ca 
20 
o 
mr 
ta 
bit 
un 
maIcn 


OPTION 2 -- 


oS Ooc2NaG @ 


amn wwNINN 
oO 


ie) 

O 

L 
s 
T M 
M TO USE T i 
s B 
E E 
N T 


HOY wM~DmamOD 


O 
s 
A 
A 
A 
U 
I 
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OPTICN 3 -- NOT DESIR 


SELECT AN CPTION: 1, 2, OR 3. 


THE sa? ha "go", “He ie pak tly 
PROM YOUR PREV 1oUS Opts 


THE FOLLOWING OPTIONS AR 
1. USE ALL OF 2 
A 


1) 
td 


E 
A 
R 
Tt 


WHS 


DISPLAYED AT 
ANGING 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLCWS:; 


{| IOL IC IR ISS IM ITF1 ITF2 ITF3 IFDFW IE IDEBUG - 
3 0 0 0 1 0 0 0 0 0 0 
ISET IDSTAE IPSD IYU INORM IREG NS NC NOB NG 


0 0 0 0 0 0 4 0 0 9 
ORDER OF SYSTEM = 4 
NUMBER CF CONTROLS = 0 
NUMBER CF CBSERVATICNS = 0 
= WUABER CF FROCESS NCISE SOURCES = 0 
‘ol 
THE SYSTEM MATRIX {"F"-MATRIX}... 
0.0 1.00000 0.0 0.0 
0.0 -0.41500 -0.01110 0.0 
§.80000 -1. 43600 -0.01980 0.0 
0.0 9.0 1.00000 0.0 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 


0.0 9-9 
8C0D+01 


ooo 
) 


0 
-0 


OPEN LOOP EIGENVALUES. ........-DET(SI-F).. 
0.0 3-6. 80767D-01: 1.22984D-01, 3. 80349D- 01: 


0.0 
0.0 2-3 
0.0 
1.900000D+00 - 


0.9 


MCDAL MEASUREMENT SCALING MATRIX...H (BAR) *T.. 
0.0 0.0 0.0 0.0 


fol, Aer ete ee ES ale BS saa te olin, tae ZEN IO ay ig cet eet ee fh Per ap eae lw iar ee Rie ee ee ee WS Wee eye TOFS. Mee oe yay Ee” at ee AO a ee ath ne, 
geri Pee yea ee PERSE IE poe ey a sa ie RO eR ee ra es Oe Oh Pa aS OTe aR seg RCT Pel on Oe ei ah PS ARE Mm 8a Sa ay aR IES RR gs Ogee ete Sets a 
BRE Geach an Peg ts, a et Ne Oe Oe ee Sale a we a i ee Se: im eee, OB ae a 


owe ie Sieh ee oh ee I Oo niger VEN SS Gg htvign Puig ite, Se i eit ene ee we! mc ea Oa: te le 0S tag ee eg ee iim Oey SR ke Se ie oe te Se Me Sie Oe Sie Sw 
SE eee LEN LE, gL ee ee pe, I Pg ee en DI I ear ty ee Bey Dee ee ee ae en NE ee ees oe pele I Gat ee eS, 


PES SEE SS SELES SS TESS SESE TE SS ESSE SES ES SESE SEES TS TS TR SEE SE ETS tS SS 


IF YOU ARE LISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 


“TYPE *Y* TO EXIT- 


(ANY OTEER INPUT TO CONTIN E) 
‘ BC ROR OE i a ok kkk 
LOADING CPTCALC.... 
EXECUTION BEGINS... 


DURING THIS SECTION OF THE PROGRAM YOU WILL 


- SELECT THE TYPE CF SYSTEM RESPONSE TO PLOT 
(CFZN LOOP, CLOSED LOOP, OR FILTER/REGULATOR) 
- PROVIDE START AND STCP f1ME FOR PLOTTING CALCULATIONS 
- SELECT THE E CF DRIVING eet ea Sone OR pepe 
- PROVIDE START AND STOP TIMES FOR TH RIVING FUNCTION (S) 
- FROVIDE DRIVING FUNCTION MAGNITUDE(S). 
CLEAR THE SCREEN TO CONTINUE 
THE F MATRIX 
0.0 1.00000 0.9 0.0 
0.0 -0.41500 ~0.01110 0.0 
$.80000 -1. 43000 ~0.01980 0.0 
0.0 0.0 1.00000 0.0 
THE G MATRIX 
0.0 
0.9 
0:9 
e es 
THE C MATRIX 
0.0 0.0 0.0 0.0 


e THE FOLLOWING PLCTTING OPTIONS ARE AVAILABLE IF THE 
FELUIRED MATRICES WERE CALCULATED IN OPTSYSX: 


1, CPEN LOOP TIME RESPONSE 
XDCT = (F}*X + {G}*UC 


+ 

2. CLOSED LOOP TIME RESPONSE 

XIOT = [{F-G*C}*X + [G}*0C, U = {C}*X 
TER 


3. Ce2TIMIZED FIL CLOSED LOOP SYSTEM RESPONSE 
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LOOP SYSTEM RESPONSE. 
Z = {H} *X 
f{*(z - H*XH}, U = [C}*XB 


AT WHAT TIME DO YOU WANT TO START 
TEE TIME RESFONSE CALCULATIONS? 


INPUT START TIME IN SECONDS. (NORMALLY 0.0) 


AT WHAT TIME DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 


INPUT STOF TIME IN SECONDS. 


RAM DIVIDES THE TIME INTERVAL YOU iAVE 
IFIED INTO UP TO 500 SMALL INTERVALS FOP 
RATION AND PLOTTING ROUTINES. IN ORDER 
OMPUTER TIME, THE NUMBER OF POINTS CAN 3E 
DUCED WITH SOME LOSS IN CURVE FIDELITY. 


HOW MANY PCINTS DO YOU WANT TO CALCULATE? 


. 


AAHGIH 


2 
500 
DOZS THE SYSTEM UTILIZE A DRIVING FUNCTION (CONT RCL INFUT)? 
(Y)ES OR (N)O 
nN 


DOES THE SYSTEM STARI WITH ALI INITIAL CONDITIONS = 0.0 ? 
(Y) ES OR (N) 0? 


THE INITIAL CONDITION FOR X( 
THE INITIAL CONDITION FOR X( 
THE INITIAL CONDITION FOR X ( 
THE INITIAL CONDITION FOR X( 


LAST OPPORTUNITY TO 
IN THE FOLLOWING AFEAS. 


ANOTHER IYPE OF SYSTEM TO PLCT 


X (3) 


FILTER OR FILTER/% EGULATCR) 
X(2) 


X (1) 


STATES AND STATE ESTIMATES CAN BE PLOTTED. 


START AND STOP TIMES 
DRIVING FUNCTIONS 
INITIAL CONDITIONS 


(OPEN, CLOSED, 
CONTINUE 


2. 
3. 
4, 
5. 
SELECT A NUMBER BETWEEN 1 AND 5. 
ALL CCNTROLS, 
0(1) 


TIME 
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24.00 0.0 -0.2804747D+00-0.8285960D-01 ).89642€11+00 
24.50 0.0 - 0.31844 75D+00-0.6822545D-01-0.5216259D+00 
25.06 0.0 - 0.347898 1D+00-0.487857 1D-01-0. 2103153D+01 


We He RF ek oe ee ee ee eK ee a kK aK RK ok RR RK eK KOK 


IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR ANT WOULD LIKE TO EXIT TOC CMS, 


¥ “TYPE 'y' TO EXIT- 
(ANY OTHER INPUT TO CONTINUE) 
PEt eS SSELESL ESS TS LTS SES SSE SS SS SSPE SSS SS SSS SSS SS SL SSS Se StS eeeLets 2 © | 


B (120) R/0 
. c (121$ k7O 
# EF (122) 5/0 ’ 
eee Your Fortran program is now being loaded ... 
eee execution will soon follow ... 
EXECUTION BEGINS...; 
e THIS PCRTICN OF THE PROGRAM PLOTS: ‘ 
- THE STATES 
- EXTERNAL CONTROL INPUTS, 
- FEELBACK CONTROL INPUTS, 
- STATE ESTIMATFS AND 
: - RECCNSTRUCTION ERRORS 
FROM THE DATA THAT YOU JUST CALCULATED. 
4 THE CAPABILITY IS ALSO AVAILABLE TO REVIEF ANY a 
: GRAPHS THAT YOU HAD PREVIOUSLY SAVED ATA 
| FILES ON YOUR DISK. 
4 CLEAR THE SCREEN TO CONTINUE. 
THE FOLLOWING OPTIONS ARE AVAILABLE: ae 
: 4. PLCI THE DATA YOU JUST CALCULATED. 
2. PLOT A CURVE THAT YOU PREVIOUSLY SAVED. 
; ENTER 1 OR 2 
| 
& 
_ YOU MAY PLOT UP TO 4 SYSTEM VARIABLES VS TIME. 
: HOW WANY VARIABLES DC YOU WISH TO 2LOT? 
. 4 
> KHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 1? 
i 1. STATE VARIABLE (IE., X1, X2, ETC) 
= 2. FFEDEACK CONTROL (If. = +c*X) 
= 32 CCNTROL INPUT (TE., U4, U2, ETC.) 
. 4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XPAT2, ETC.) 
2 5; STATE RECONSTRUCTION ERROR (IE-., X1-XhaTi, 
. X2-XPAT2, ETC) 
. ENTER 1,2,2,4 OR 5 
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WHAT IS THE SUBSCRIPT OF THE STATE VARIAPLE THAT 


ry YCU WANT TO PLOT AS THE NUMBER 1 CURVE VS TIME? 5 
1 § 
WHAT IS THE CURVE LABEL FOR THIS VARIABLE? - 
Ec NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
F) => BETA 
F) => PHI 
0) => THETA 
2 state yl 
wld 
WHICH TYPE OF VARIABIE DO YOU WISH TO PLOT AS CJEVE NUMBER 2? 
J. STATE VARIABLE (15, cETO) 
2. FEEDEACK CONTECL (if. Goce 
32 CCNTROL INPUT 025 ETC. 
4. STATE ESTIMATE tates, (1E., HAI1, XHAT2, ETC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-XMaTI1, 
K2-XHAT2, bee 
P ENTER 1,2,3,4 OR 5 
| 
a KHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT al 
‘ YOU WANT TC PLOT AS THE NUMBER 2 CURVE VS TIME? 
@: 2 
: WHAT IS THE CURVE LABEL FCR THIS VARIABLE? 
d NOTE: 1. 40 CHARACTERS MAX LENGTH 
- GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IF. (A) => ALPEA 
E} => BETA 
F) => PHI 
Q) => THETA 
j state y2 
WHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CJRVE NUMBER 3? 
1. STATE VARIABLE (IE. X1, X2, ETC) 
2. FEEDEACK CONTEOL (TE. =c*X) 
3° CCNTROL INPUT (IE., U4, U2, ETC.) 
G2 STATE ESTIMATE (OBSERVER) (1E., XHAT1, XHAT2, ETC.) 
- S. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC) 
S Z ENTER 1,2,3,4 OR 5 
2 1 
p 
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HARACTERS MAX LENGTH 

K SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
OSED IN PAKENTHESES. 

> ALPHA 


> BE 
> PH 
> T 


Hta 


HI 
HETA 


e 
—_—— 
or4 tri ee 
Peet 
eu Ut fll 


state y3 


WHICH TYPE OF VARIABIE DO YOU WISH TO PLOT AS CJRVE NUMBER 4? 


1. STATE VARIABLE (IE., X1, X2, ETC) 
2. FEEDEACK CONTEC f. “c*X) 
32 CCNTEOL INPUT i vi, gos ETC. ) 
4. STATE ESTIMATE oBserRvé ER) (LE, HAT1, XHAT2, ETC.) 
5:2 STATE RECCNSTRUCTION ERROR (IE:, X1-xhaT1, 
X2-XHAT2, ETC) 
: ENTZR 1,2,3,4 OR 5 
? 
“WHAT IS THE SUBSCRIPT OF THE STATE VARLABLE THAT 
‘ YOU WANT TC PLOT AS THE NUMBER 4 CURVE VS TIME? : 
4 ~ 
WHAT IS THE CURVE LABEL FOR THIS VARIABLE? : 
NCTE: J. 40 CHAEACTERS MAX LENGT 
3. GREEK SENBORS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IF. (A) => ALPHA 
P) => BETA 
F) => PHI 
C) => THETA 
state y4 
YCU MAY USE UP TO 3 HEADINGS. 
HOF MANY HEADINGS DO YOU DESIRE ON THIS GHAPH? 
“ 0, 1, 2 OR 3 
3 


WHAT IS THE DESIRED HEADING NUMBER 1? 


NCTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREE YMBOLS WILL BE PR ED FOR ANY 
LETT ESES. 
IE. => AL 


R 

S) INT 
S ENCLOSED IN PARENTH 
| PHA 


=> THETA 
open loog systen 
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WHAT IS THE DESIRED HEADING NUMBER 2? 


{ NOTE: 1. 40 CHARACTERS MAX LENGTH ’ 
2. GREEK SYMBOLS WILL Bz PRINTED FCE ANY 
LETTERS ENCLOSED IN PARENTHESES. 
IEF. (A) => ALPHA 
B) => BETA 
F) => PHI 
: Qo) => THETA 
| example 2 ’ 


WHAT IS THE DESIRED HEADING NUMBER 3? 


NOTE: a 40 CHARACTERS MAX LENSTH 


LETTEKS ts ane IN PARENTHESES. 
: IE. (A aN PHA ’ 
B) => BETA 
F) => PH 


me 
Q) => THETA 


modern centrol theo 
>> USING A PRE- ALLOG TED DATASET FOR UNIT FT17F001. 
TEES FOLLOWING OPTIONS ARE AVAILABLE. ; 
: BEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. % 
2 REPLCT FREVIOUSLY SAVED GRAPH DATA. 
3. EDIT THE CURRENT GRAPH. 
- FLOT REVISED GRAPH ON THE TEK618. 
5 pole ae eae MAKE METAFILE OF THE CURVES. 
REVIOUSLY SAVED. 
; SELECT A NUMBER BETWEEN 1 AND 5S. ’ 
3 
THE GRAPH EDIT MEND 
1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 
- 2- DELETE CURVE FROM CURRENT PLOT. ” 
J 3. EDIT CURVE goer eats ' 
4. EDIT PASE HEADING(S). 
5. CHANGE THE Y-AXIS SC CALE. 
6. CHANGE THE TIME AXIS SCALE. 
7. CHANGE PLOT SIZE. (Be Faure Is 8.5 X 6.0) 
8. CHANGE THE LETTERING HEIGHT. 
9. CHANGE POSITION OF THE LEGEND. 
» 1). EDITING COMPLETE. , 
SELECT A NUMBER EETWEEN 1 AND 10. 
i 
HOW MANY INCHES IN THE X DIRECTION 
LEFT OR RIGHT), DO YOU WANT TO MOVE 
> GOVE THE LEGEND BOX FROM ITS PRESENT PISITION ’ 
NOTE: 1. DEFAULT FLOT SIZE IS 3.5 X 6.0 
2. LEFT IS NEGATIVE 
a 3. RIGHT IS POSITIVE 
-4 
2 HCW MANY INCHES IN THE Y DIRECTION ’ 
; sr OF DOWN) D YOU FANT TO MOVE 
MOVE THE LEGEND BCX FFOM ITS PRESENT P)SITION 
37 
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8.5 X 6.0 


THE GRAPH EDIT MENU 


ooo 


(8 CHARACTE 
S 
6 
R 
E 


e . ry C 
Ww Yn Y [oe] 
. fj &y in hy] 
& a] is) oO AQ 
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Ay H H H 
toma] te Uba A 
ica] 2p) 9) fag WwW 
=z is) iQ Hes [a9 
&] oO «ix (je) 
fs} fae) m~ a] n= EY 
fee) a fo) fay Yn 
D me fyett 
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oO n un «t fx ot 
fx e fq ry 4 e a Te C4 
7 oO & n ios zp) Eu«e (=) 
H rc fq << fx ee Om 
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THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAFILE OF PREVIOUSLY SAVED IURVE. 
2. QUIT. 


ENTER 1 OR 2 


WHAT FILE NAME IS THE DATA STCKED UNDER? 
openloop 


THE CURVE DATA I 


1 S_BE LO 
>> USING A PRE~ALLOC D 


NG A 
IED DA 
THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAFILE OCF PREVIOUSLY SAVED CURVE. 
2. QUIT. 


ENTER 1 OR 2 
- 13197 VECTO 


~- 1 KS 
TWARE PRODUCT OF 
E REFERENCES; 5 


GENERATED 
ISsCO, SA 
EADS;'0 W 


1 PLOT FRAMES 
5GO, CA. 
R S. 


WON bito-v 
Pol eslestades1o) 
mean Mm 


He ee he he ee ae oe oe oe eae ee te ee a oe Re ae oe a a Re a ak ak eo ooo kk ok ak ee 


DO_ YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
TEE DISSPLA ETA oR Ay YOU JUST CREATED? 
r\ 


eRe he oe eo a ke oe ae eo ee ocak ke ie ka a ca ee eo eo ea te oe co ake eke teak ea koi oe eo ae otek oe oe fe 


R/O; K/W BY MVS 


MVS COPY 091 NOHOLD 


ODMAN 


We ee he ie he eee he eee eae Re ae a oe I ea a 2 he Re ot a ke ee he ie eae Ra a Re ak ke ok ke a ke ee ke oe 


YOUR GERAPH(S) CAN BE PICKED UP AT THe COMPUTER CENTER. 
TEE GRAPH (S) WILL BE ALDRESSED TO "POP (JSERK ID)". 
2 GG Rk lok Rito gto ok tik fo ago do a Rk do dk tok soko ek 


es He Ree ee ee ke ee ee ee Ee Re ae a ea a te a a ee eo a ok ak tee trate ok 


DO YOU FANT TO 


1. FUN OPTSYSX AGAIN 
ze RUN THE FLOT PROGRAM USING THE SAME MATRICES? 
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eras 
eer? 


Sek ghcE tee Be der ehe 4 See ea We poeta ele 


( 


eGre ANOTHER TYPE OF SYSTEM (OPEN/CLOSED) ) 


ENTER 1, 2 OR 3 
3 RG ido kok loki dk tok kkk ik fk kk ok kek oto kok ok tok a tok dk ok ok 
3 


KH Re RH RR Re Re I a ee KK KK a eK eK KK 


HAVE A GOOD DAY!! 


Tce SSS S ESS SC SPECS SSS SESS SS SS ESS SSS SSS SESS SSS SPS CSTR SET SS SSS tS SS 
Rs T=11.05/17.26 20:08:15 

record off 

END RECCEDING OF TERMINAL SESSION 


The grazhical output generated by this exanple follows 
as figure 3.1. 
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sett 


OPEN LOOP SYSTEM 
EXAMPLE 2 
MODERN CONTROL THEORY 


TO- 


TA ULVLS 
(ete eae ee er a ae hie an ee alte ae a a eee eel 
0sa°0 $20°0 000°0 $z0'0- 00°0- §L0°0- 00t'0- 

eA TLVILS 

‘e@ 

Ce ee ee ee 
o's o'r o's oz ot oo Ot- oz- = O'E- 

EX GLVLS 
0°02 o'st oot a9 0° o'9- aOt- 

PAGLVLS 


4] 


SORE tet ARB ts el re OS is Rae onan teed te Nana! iat aie Coss ec ere) oe gm 
BT te ee aR ee ee Pe weg eT ek ete wie fore ey ere et kennels 


TIME — SEC 


Open-loop Time Response 


Figure 3.1 


EB. CICSED-LOOP SYSTEMR TIME RESPONSE 


The foliowing closed-loop system exanple was taken from = 
{Ref. 6, pp 5.8 - 5.19]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case Letters or numbers 
below each "2", a 


‘wo 


ING OF (eas SESSION 
02 20:19:44 


1. OPTSYSX FCRIRAN SYSTEM ENA RES) 
2. OPTCALC FCRTRAN CALCULATE TIME eee 
3. OPTPLOT FCRIRAN DISSPLA PLOTTING ROUTINE) 
EACH PRCGRAM PASSES INFORMATION TO THE NEXT PROGRAM 
THROUGH A DATA FILE WRITTEN TO THE USERS DISK. IN THIS 
CASE, THESE FILES ARE “OPTMAT DATA" AND “OPTPLOT DATA", 
THE SIZF OF THESF FILES VARY WITH -THE SYSTEM ORDER, AND 
CAN USE ABOUT 20% OF THE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 
- TYPE "E" TO EXIT, ANY CTHER ENTRY TO CONTINUE - - 


Be ee ee ee ee a he he ee Re a oe ie a eo ee ae oe 2 ee a he a oe fe fe ek tate ea oe ok ake a ako kK ae ak 


YOU HATE A DATA peeve NAMED Sveeaee DATA! oF YOUR A 


DISK THAT Ha VIOUSLY GENERATED BY TH ee e 
EFROGRAM AND RTAINS THE F, ee H, GAMMA, A AND - 
MATRICES FROM THAT RUN. nove 
IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 

THIS _KUN CPTSYS PROGRAM WILL READ IN THE 

CESIRED DATA AT THE APPROPRIATE TIME 


f 


IF YOU TYPE (Y) ES. 


ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 


PLS ELE LE SES ET TES ST SS SEES STS SS TSS TS TS PSSST TSS TS TS SS SRE SS SS tt SS 
Y 


Se He Ke ee oh eee ee ee ke ee ae oe oe oe ok Re eo te RE a ek Kk & 


CC YCU WANT THE NUMERICAL CUTPUT FROM OPTSYSX TO GO 
TC YOUR TERMINAL (bE DS OR TO A D(ISK) FILE? 


Me Ke He te ee Ree He eK hea RR eo a eo ee eK a ok ee kK 
s 


PES SRE SESE SESE ES ETS SSE SEL TS BESS ESSE SCE SS STS SCS TESS SSE ST SETS tS SS 


OUTPUT WILL COME TC YCUR TERMINAL SCROUCH. 
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PEL SES ELE SES TESS ETE SESE SE SRS SESE TES ESE TET ESS TRC SS ESS SES ET ETE SS 


EXEC EGINS, 


Lente fer lee) 
ot 
mw 
Pad 
oe) 
| 
wi 
< 
a) 


Attn ht 
Or 
1] al) 
Ke 
OuiH 
Y 2 
ty2H 


6) 
MEROUS CO 
MS CONTROL 


= {F}*X + [G}*U + {GAM} *(W+HO) 
MEASUREMENT EQUATION-- 
{H}*X + {D}*U + V 
REGULATOR PERFORMANCE INDEX-- 
1/2 * INTEGRAL (Y *{A}*Y + U *{B} *U) DT 
STATE FEEDBACK GAIN DEFINITION-- 

U = - {C} *X 
DC YOU WISH TO CONTINUE? TYPE "YES" CR "won", 

--DATA ENTRY-- 

SPECIFICALLY DE 
ie TO USERS: 


uP a see Gu AND "GAMMA" 
L FROM SEPARATE DATA F 


ura Ud gu A 
DEFINED’ WI 
A 


N 
u 


Qa 
tt 


HQ 


SI 
KA 


EXPLICI 
SE, TH 

THE PROGRAM’ LIS 

THE EXAMPLES CONTA 
I 


E 
T 
I 
DO YOU WISH TO CONTIN 


I 
N 
U 


rO YOU WISH TO INPUT THE "F", "G", AND "GAMMA" 
MATRICES FECM SUBROUTINE “SETUP" IAW THE 
METHOD DESCRIBED ON THE PREVIOUS SCREEN 
TYPE "YES" OR "NO". 


GENERAL CPTSYSX CPTIONS: 


oh) 


CFTION 1 ~- SYSTEM Bees WITHCUT 
CPEN-LOOP EIGENSYSTEM CALCULATIONS. 
CPTION 2 -- SYSTEM ANALYSIS WiTH OPEN-LOOP 
FIGENS YSTEM CALCULATIONS. 
CPTION 3 -- CPEN-LOOP EIGENSYSTEM FOUND 
ND PROGRAM TERMINATES. 
{UF"-MATREIX ENTRY FOLLOWS IMMEDIATELY. } 
CETION 4 == MCDAL DISTRIBUTION MATRICES COMFUTED 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
CR STEADY-STATE ANALYSIS. 
SELECT AN OPTION: 1,2,3, OF 4. 


Bile Tw OED, Sea ne tue Ee Tg SORES, ae ye tg Ee tee Ee Ree ye Tee OR SAS an Ld Sena Le A Sag Aa SQL eae ON SRE OM SSS Re Se oS 
ee ee eS ae, Tet Se Fel isi Fe ie a a em ne eS mo ce ig ae ig Ma Pe a, Biel Me Sgt ae tg ta Se a 


DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
TYPE "YES" CR "NO", 


n 
CETSYSK LQOR/CLASSICAL OPTIONS : 
CEFTION 1 -- CPTIMAL FILTER AND/OR REGULATOR 
SYNTHESIS WITH NO EXTERNAL "C" OR “K" 
MATRIX INPUT. 
CETION 2 -- CPTIMAL FILTER AND/OK REGULATOR 
SYNTHESIS WITH EXTERNAL "C" 
MATRIX INPUT. 
CETION 3 -- OPTIMAL FILTER EAS REGULATOR 
SYNTHESIS WITH EXTERNAL "K" 
MATRIX INPUT. 
CETION 4 -- CPTIMAL FILTE aeeoe REGO 
SYNTHESIS WITH EXTERNAL "C" AND “K" 
MATRIX INPUT. 
S SELECT AN OPTION: 1, 2, 3, OR 4. 
4 
DO YOU WISH TO CETERMINE THE STEADY-STATE RESPONSE 
FOR A CONSTANT DISTURBANCE? 
TYPE "YES" CR "NO", 
DO YOU WISH TO LETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 
TYPE "YES" CR "NO", 
n 


CEEN-LCOP TRANSFER FUNCTION OPTIONS: 


CPTION 1 -- NO OPEN-LOOP TRANSFER FUNCTIONS CONPUTED. 
CFTION 2 -- FOLES, RESIDUES, AND ZERCS COMPUTEL. 
CETION 3 += CNLY POLES AND ZEROS COMPUTED. 

CFTION 4 -= CNLY PCLES AND RESIDUES COMPUTED. 


SELECT AN OPTION: 1, 2, 3, OR 4. 


NCISE TRANSFER FUNCTION OPTIONS: 


CETION 1 -- NO NOISE TRANSFER FUNCTICNS COMPUTED. 
CFIICN 2 -- EFOLES, RESIDUES, AND ZERCS 2OMPUTETD. 
CETION 3 -- ONLY POLES AND ZEROS COMPUTED. 

CETION 4 -- CNLY POLES AND RESIDUES COMPUTED. 


SELECT AN OFTION: 1, 2, 3, OR 4. 
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= 


ow 


COMPENSATOR TRANSFER FUNCTION OPTIONS: 


OFTION 1-- NO CCMP. TRANSFER FUNCTIONS COMFUTED. 
CETION 2 -- FOLES, RESIDUES, AND ZEROS COMPUTED. 
CPTION 3 -- CNLY POLES AND ZEROS COMPUTED. 

CETION 4 == CNLY POLES AND KESIDUES COMPUTED. 


{NOTE: A COMPENSATOR TRANSFER FUNCTION CAN BE 
COMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
AND/OR INPUT. 3 


SELECT AN OPTION: 1, 2, 3, OR 4. 


WILL A FEED-FORWARD DISTEIBUTION MATRIX 
{"D" - MATRIX} EE INPUT ? 


TYPE "YES" OR "NO", 


mraH 

ctttioor 

Mie S@rIDH 
sire ae DD 
wn Ode 


aru wHO 
_—- Wom 
D HRH 
~~ oO YE 


I 
0 
RENT MARKOV CRITERIA, 


im ow 


E LEFAULT VALUE, TYPE "0" {(ZER 0} 


DO YOU DESIRE TO SYNTHESIZE A STABLE FILTER {CR KEGULATOR} BY 
DESTAEBILIZING TEE ORIGINAL SYSTEM? 


{NOTE: WORKS FOR FILTER OR REGULATOR BUI NOT FOE EOTH 
IN THE SAME RUN.} 


TYPE "YES" CR "NO", 
[Tc YCU DESIRE TT’ PRINT THE EULER-LAGRANGE EIGENSYS‘TEM 
FRICK TO DECCHPCSITION {FOR CHECKING THE PROGRAM}? 
TYPE "YES" CR "NO", 


PCWER SPECTRAL DENSITY {PSD} OFTION 1 : 
OFTICN 1 -- COMPUTE THE PSD OF “HE OUTPUTS AND/CR THE 
CONTROLS OF THE CONTROLLED SYSTEM WEEN FCFCED BY 
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tte sen ee si ee ee Oe ete See Bie a mie Sie ee Rie rN cS ee ant diigt t ey ®: 


PROCESS AND MEASUREMENT NOISE. fecae: BOTH A 
REGULAICR AND A FILTER MUST BE RESIDENT IN THE 
PROGRAM TO USE THIS OPTION.} 
‘t OPTION 2 -- SAME AS OPTION 1 ABOVE BUT ONLY PRINT THE ' 
RESIDUES OF BACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 
OPTION 3 —- NOT DESIRED 


SELECT AN GCETION: 1, 2, OR 3. 


Gaed 


DC YOU DESIRE REGULATOR SYNTHESIS ONLY? 
TYPE "YES" OR "NO", 


THE "Fe, ge, wy, uGaMt, MAM AND "EM WATRICES 
FROM YOUR PREV IOUS OPTSYS RUN WERE SAVED. 


e THE FOLLOWING OFTIONS AKE AVAILABLE: ' 
1. USE ALL OF THE SAME MATRICES ASAIN. 
- USE SELECTED MATRICES AGAIN. 
3. INFUT ALL NEW MATRICES. 


ENTER 1, 2, OF 3. 


NOTE: Rea SAVED MATRIX WILL BE REDISPLAYED AT Seek, 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
4 INDIVIDUAL MATRIX ELEMENTS. : 
2 
r CO YOU WISH TO SAVE THE "FY"-MATRIX FROM THE LAST ee 
ta KUN TO BE USED IN THIS RUN? 
NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE FRCPER INPUT SEQUENCE INTERVAL 
AND YOU VILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
TYPE "YES" OR "No". 
ed 
Y 
[CC YOU WISH TO SAVE THE "A"-MATRIX FROM THE LAST 
FUN TO BE USED IN THIS RUN? 
NCTE: THE MATRIX WILL BE REDISPLAYED AT 
e THE PRCPER INPUT SEC pee INTERVAL 
AND YOU WILL HAV PTION OF CHANGING 
INDIVIDUAL MATRIX ELEM ENTS. 
TYPE "YES" OR "NO", 
Qn 
CO YOU WiSH TO SAVE THE "B"~MATRIX FROM THE LAST 
FUN TO BE USED IN THiS RUN? 


46 


Regrt iS ele SENG Seite Chats vet taetinn Supa “Rete Sah ky RTS a ae emia SU. SR ORR RS gree, ea ROR LOT Saye agua te tS a oa 

Poe ee” a et Bee Sree. Ra! ea Jee. aoe wey seo alae ete. Sens er ttin OSE Am ya lh eo, ee MiBb arte y leet | GR Goer ei 8 een gS Fes PE here Cen, I Sober, Ge et he 

ee oe gee ea taht el acme ale We ee Ute ee Lene er cby tty Medge Renee RON ge Sn wate. Co or Regattas eis oh eter a Beton cea. Seley ORR Sm a Rye eran e See ime ge ae A ok ee ee ae, A ae 

ge gi ee a a ee Oe we eee me eee, ene te ee te ee, ONS He. ie Nn ieee era Mea cee Sete aE we Mae Page ae eat ge) gt Seg! Sask i ap eat ag aS Se bata oe Se a 
a eee nee Oe 


aS ak iat tp ee Sa Gel we eo a 


NCTE: “HE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE TRE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
TYPE "YES" OR "No" 
n 
ENTER THE # OF CONTROLS {NC} OF THE CONTROL SYSTEM MODEL 
,  ("GN-MATRIX}. 
1 
ENTER THE f OF MEASUREMENTS OR OBSERVATIONS {NO} OF THE 
("H"-MATRIX}. 
? 
rf 
ENTER THE # OF FROCESS NCISE SOURCES {NG} OF THE 
‘ iNGAM MA MATRIX}. 
0 
FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLIC¥S: 
IOL I¢ IR ISS IM ITF1 ITF2 ITF3 IFDFW IE IDEBUG 
9 0 0 Oo 0 0 0 0 0 0 0 
ISET IDSTAB IPSD IYU INORM IREG NS NC NOB NG 
0 0 0 0 0 1 4 1 4 0 
ORDER OF SYSTEM = 4 
NUMBER CF CONTROLS = 1 
NUMBER CF CBSERVATICNS = 4 


NUMEER CF FROCESS NCISE SOURCES = 0 


THE SYSTEM MATRIX {"F"-MATRIX}... 
Q.9 1.00000 0.0 0.90 
0.0 -0.41500 -0.01110 ° 
$.80000 -1.43000 -0.01980 0.0 
0.0 0.0 1.00000 0.0 


DO YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


TYPE "YES" CK "NO". 
D 
CEEN LOOP DYNAMICS MATRIX. .cccsccccesccaaeFes 

0.0 0.10006+01 0.0 0.0 

0:0 -0:415(D+#00 -01110D-01 9: 

0.98COD+01 -0.14305+01 -0.1980D-01 0. 

050 0.0 0:1000D+01 0. 

ENTER THE MEASUREMENT SCALING MATRIX {"H"-MATRIX}. 
DIMENSION = #. OBSERVATIONS {NO} X # STATES {NS} 

THE ELEMENT H( 1, 1)= 
j 
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sitighie te San? se 
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ELEMENT 


FLEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


FLEMENT 


ELEMENT 


ELEMENT 


H( 1, 2) 
H( 1, 3) 
H( 1, 4) 
H( 2, 1) 
H({ 2, 2) 
H( 2, 3) 
H( 2, 4) 
H( 3, 1) 
H( 3, 2) 
H( 3, 3) 
H( 3, 4) 
H( 4, 1) 
H( 4, 2) 
4H( 4, 3) 


H( 4, 4) 


THE MEASUREMENT SCALING MATRIX {"H"-MAT RIX}... 


0000 


0.0 
00 0:0 
1.00009 


TO YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO". 


MEASUREMENT SCALING MATRIX... ces ceceoces eH. 


THE 
? 
0 
THE 
? 
0 
, THE 
: ? 
y 0 
: TEE 
: % 
] 0 
THE 
? 
1 
THE 
? 
0 _ 
7 ? THE 
. 0 
THE 
? 
0 
THRE 
? 
0 ™ 
» “ THE 
4 
is THE 
) 
THE 
~ ? 
0 
. THE 
0 
THE 
2 
: THE 
= E 
ee 
1 
029 
» 0.0 
0.0 
ry nN 
eh 9.1 
‘ 0.0 
v 0.0 
® 0.0 
S 


0core401 


000L+01 


0.9 
0.0 
Q. 1000 
0.0 
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0.0 
0.9 


D+01 0.0 
0.1000D+01 
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OV OV OV OV OV OF OV OV OV OV OV 


Od OFF OFF OV Ov 


e 
tS 
(6) 


CUTEOUT 


9.0 


ENTER Ud MEASU 
DIMENSION = # OESERVAT 
THE ELEMENT A( 1, 1)= 


THE 


THE 


lolelole) 


OOOO 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


THE OUTP 


ELEMENT A( 1, 
ELEMENT A( 1, 
ELEMENT A( 1, 
ELEMENT A( 2Z, 
ELEMENT A( 2, 
ELEMENT A( Z, 
ELEMENT A( 2, 
ELEMENT A( 3, 
ELEMENT A( 3, 
ELEMENT A( 3, 
ELEMENT A( 2, 
FLEMENT A( 4, 
ELEMENT A( 4, 
ELEMENT A( 4, 


ELEMENT A( 4, 


MATRIX}. 
S {NO 


OUTPUT MEASUREMENT COST MATRIX {"A"-MATR IX}... 


0.90 
0.9 
Q.0 
0.0 


TYPE "YES" CP 


0.9 


0.0 
QO. 
0 


oo 


WHOM, 


0.0 


4g 


9.0 
-0 


0 


Ooo 


COST MATR oY Xo vccee 060 0's 6.0% 000 6 6 eiee sete Ase 


0.0 


«25000 


9.0 0.0 0.0 0.0 
0:20 0.0 0.0 .0 
0.C 0.0 0.0 0.2500D+00 
i ENTER THE CONTRCL DISTRIBUTION MATRIX {"G"-MATRIX}. 
DIMENSION = # STATES {NS} X # CCNTROLS {NC} 
THE ELEMENT G( 1, 1)= 
0 
5 THE ELEMENT G( 2, 1)= 
l 6.27 
_ THE ELEMENT G( 3, 1)= 
$.8 
i THE ELEMENT G( 4, 1)= 
. 0 
b 
THE CONTROL DISTRIBUTION MATRIX {"G"-MATRIX}... 
0.0 
6.27000 
S- 80000 
0.0 
) 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO", 
n ; 
r ENTER THE CCNTRCL COST WEIGHTING MAT TEI "BM MAT 7a aa 
" DIMENSION = # CONTROIS {NC} X by ROLS ine 
: THE ELEMENT B( 1, 1)= 
131.3 
TEE CONTROL COST MATRIX. .ccccccee Boos 
6 131.30000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYP "YES" CR "NON, 
nN 
» 
THE CCNTROL DISTRIBUTION MATRIX... c...2e00eGee 
0.0 
9262705401 
9. 98C0D+01 


TEE CONTROL COST MATRIX ccc ecw cc ee ccc snes ce Bee 
0. 13130+03 


EIGENSYSTEM OF OFTIMAL KEGULATOR......0.- 
C-LOCP CPTIMAL FEG. E-VALUES...DET(SI-F+G*C).. 
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—1.23385D+0C, 5.54546D-01:-4.198350-01, 1.13532D+ 00: 


C-LCCP RIGHT EIGENVECTOR MATRIX... ...2200eMeeee 
-1.019344D-01 08717D-02 -8.155484D-02 1.1222 64D-Q1 
1.12969 1D-01 as 02 -9.317336D-02 -1.397) 74D-01 
1.0000 90D+00 “3: 1.000000D+00 0.0 
-6.742684D-01 -3. $30447D- 01 -2.665351D-01 -7.7484S59D-01 
CCNTROL EIGENVECTICR MATRIX... ccc ce eee ee CHM. 
-5.4€4374D-03 2.109409D-02 2.7139259-02 -1.676334D-02 
C-LCCP OPT. REG. LEFT E-VECTOR MATRIX. .M-INV.. 
-3. 764753D+00 2.578703Dt0Q0 -3.562309D-01 -1.010220D+00 
-3.421605D+01 -9.486653D+00 -4.604269D+00 -3.245261D+00 
3.764753D+00 -2.578703D+90 1.356231D+#00 1.010220D+00 
1.526581D+#01 2.419€63D+#00 1.609198D+00 4.841548N-01 


THE CFETIMAL FEEDBACK GAIN CONTROL MATRIX. ..C=BINV*GT*S... 
~8.54920-01 -3.2475D-01 -8.5345D-02 -4.3635D-02 


DO YOU WISH TO OBTAIN A TIME RES 
OF THE SYSTEM YOU ARE EVALUATI 

OR N) 
NOTE: YOU M 
(TEK 6 
THE ar gozSze™) 
( Cc 


y 


He Re eR a ee eee he eae he ee ee ee ee ake he he a he ee ee ok ke eta ie ke oe ti ok kak ak ok 


IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR ANT WOULD LIKE TO EXIT TO CMS, 


“TYPE "y' TO EXIT- 
OTHER INPUT TO CONTINUE 
% *& 


eeneseenavense nes KEKE KE RKEEKEKEK becdeeehieueaeens cea 
LOADING OPTCALC.... 
EXECUTICN BEGINS... 
CURING THIS SECTION OF THE FROGRAY YCU WILI 
- SELECT ane a6 CF SYSTEM RESPONSE TO PLOT 
OPEN LOOP OSED LOOP, OR FILTER/REGULATOR 
ROVIDE STA RE ANE STOP f1ME FOR PLOTTING CALCULATIONS 


a ee ati aie SWS iat Me cae Fa op J a SSD ee Math TN oat RS Bias Reg 1%, aE tS SSNS, Caimi cacy MSE Se yee tee aes Miah ah gs ase ose gs 
ap i aig Stee OS ERS Ree Reais op, cee tegen Ma fF ere Ee ee te eta Se hele Sena a Weer hehe, aren Sur tae MES oe 
ace cant tet On yee Oe Me ght ag a ee eR lege) Og eee Ea ater im etre aig Me ILS ge cai rg a Oe moe tn eee, Sta ery ae at Sete ote inks Oe ee Ne oui im WI ee ETN ie en ete Ake 


= SELECT THE TYPE CF DRIVING FUNCTION (S) (STE2 OF RAM?) 
- EROVIDZ START AND STOP TIMES FOR kIVING FUNCTION(S) 
: - FROVIDE DRIVING FONCTION MAGNITUDES) - 
¢ CLEAR THE SCREEN TO CONTINUE ' 
; THE F MATRIX 
I 0.0 1.00000 0.0 0.0 
0.9 -0241500 -0201110 0:0 
ke §.30000 -1. 43000 ~0:01980 0:0 
0:9 0:0 1.00000 0:0 
y 
THE G MATRIX 
: 0.0 
¢ 6.27000 
¢ §. 80000 
0:90 
THE C MATRIX 
-0.85462 -0.32475 -9.08534 -0.04364 
THE H MATRIX 
1.00000 0.0 0.0 0.9 
0:0 1.00000 020 920 
0:9 0:0 1.00000 0:0 
0:9 0:0 0:0 1.00009 
THE K MATRIX 
0.0 0.0 0.0 0.0 
0:9 0.0 0.0 0:0 
0:9 0:0 0:0 3:9 
0:0 0:0 0:0 020 


THE FOLLOPING PLOTTING OPTIONS ARE AVAILABLE IF TH 
RKEVJIRED MATRICES WERE CALCULATED IN CPTSYSX: 


1.  C2EN LOOP TIME RESPONSE 


ts 


KDOT = {F}¥*X + [G}*UC 
2. CLCSED LOOP TINE RESPONSE 
XDOT = {F-G*C}#X + {G}*#0C, U = {C}*X 
3. OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 
XDOT = _{EL*X + (6) OC, {hh 
XHDOT = {F}*XH +#{G}U + {K *(2 - H*XH} 
4. FILTER + REGULATOR CLOSED LOOP SYSTEM PESPINSE. 
ALOT =_(EtGeC} 4K + (6) #0C Z = ae 
XHDOT = {F}*XH +{c})U + {K{*{Z - H*XH}, J = {C}*XH 
es SEIPCT “1,: 2). 3-0R- 4. 
2 
THE AUGMENTED F MATRIX (F+G*C) 
0.9 1.00000 0.0 220 
-£2 3€034 -2.45120 -0:54621 -3.27353 
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1.42180 4.61257 - 85618 -0.42763 
0.0 0.0 - 00000 0.0 

AT WHAT TIFKE DO YOU WANT TO START 

TEE TIME RESPONSE CAICULATIONS? 


INPUT START TIME IN SECONDS. (NCRMALLY 0.0) 


AT WHAT TIME DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 


INPUT STCP TIME IN SECONDS. 
25 


4A 
T S 
EATION AND PLOTTING ROUTINES. IN OR 
OMPUTER TIME, THE NUMBER OF POINTS C 
DUCED WITH SOME LOSS IN CURVE FIDELi 
, MANY PCINTS DO YOU WANT TO CALCULATE? 

500 

DOES THE SYSTEM UTILIZE A DRIVING FUNCTION (CONI ROL INPUT)? 
(Y)ES OR (N)O 


n 
DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 
: (Y)ES OR (N)O? 
: WHAT IS THE INITIAL CONDITION FOR xX( 1) ? 
0502 
: WHAT IS THE INITIAL CONDITION FOR X( 2) ? 
6 
4 WHAT IS THE INITIAL CONDITION FOR X( 3) ? 
) 
? WHAT IS THE INITIAL CONDITION FOR X( 4) ? 
0 
MES LGnIQUE HPP AQ°EQEIBITAE BPs. 
1, SELECT ANOTHER TYPE OF SYSTEM TC Pio 
(OPEN, CLOSED, FILTER OR FILTEE s EGULATOR) 
2. START AND STOP TIMES 


3. DRIVING FUNCTIONS 
4, INITIAL CONDITIONS 
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CONTINUE 


5. 


° 
wn 
(=) 
= 
<< 
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= 
4 
wy 
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& 
WH 
[oe] 
0G 
fq 
mw 
= 
2p 
= 
x 
& 
O 
fQ 
Hv 
Ww 
YN 


PLOTTEL. 


CAN BI 


STATES AND STATE ESTIMATES 


ALL CCNTROLS, 
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elolojelelelelelelelelelelelalsilelelolelelelelelelelslelelololelelelolelololelelelelsiolelelelelele) 

ANOMAONOMOMOWMOMO NOMNOMNOMNDMNANDNIONONOMDNONONANOMOMONOMGE 
eoee74rtrteereteteesvsreesteoaseeeseeeseeesseeeseeeeseeereeteeeereeeeeeervseesee 6 
DOPE NANAMNM ST FUMNDOVE MH DONHNOO LTH ANMMT FNMNOVEEODANGOTEANMNM SAWN 
FTF EHH EK KM NNNNANNNANAGN 


Me He He HH HH HHH eK ee ee eke he eo eo eK RE OR Re 


IF you ARE LISSATISFIED WITH TH 
THUS FAR AND WOULD LIKE TO EXIT 


-TYPE 'y* TO EXIT- 
(ANY OTHER INPUT TO CONTINUE) 
$k a ok kk do ok do ok ok 


0) R 
| Ye 
2 


Ri 


x 
2 
2 
2 


Sok Rk 
+ 
Cc 
E 
«ee Your Fortran program is now being loaded ... : 
~e. Execution will soon follow ... | 
EXECUTION BEGINS... 


THIS OF TRE PROGRAM PLOTS: | 
TATES 

NAL CONTROL INPUTS, 

ACK CONTROL io 
FSTIMAZTES AND 

STRUCTION ERROR 

A THAT YOU JUST CALCULATED. 
T 


Y IS ALSO AVAILABLE TO REVIEW ANY 
YOU HAD PREVIOUSLY SAVED AS DATA 


SPSS cami alas a oe 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. PLCI THE DATA YOU JUST CALCULAI ED. 
2- PLOT A CURVE THAT YOU PREVIOUSLY 


SAVED. . 
ENTER 1 OR 2 a 


YOU MAY PLOT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DC YOU WISH TO PLOT? 


4 


VHICH TYPE OF VARIABIE DO YOU WISH TO PLOT AS CURVE NUMBER 1? 


1. STATE VARIABLE (IE., X1, X2, ETC) 

2. FEFDBACK CONTRCL (If, f = X) 

3. CCNIROL INPUT (IE, ile bo. ETC.) 

4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAT2, FTC.) 

5. STATE RECONSTRUCTION ERROR (IE., X1-XHaA™1 

K2-XHAT2, ETC.) 
: ENTER 1,2,3,4 O8 5 

WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TO PLOT AS THE NUMBER 1 CURVE VS TIME? 
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EO, Or, Wee Snear Wn Seay Sameer 


WHAT IS THE CURVE LABEL FORK THIS VARIABLE? 


NOTE: 1. 40 mec eeaaere EAD LENGT 
2. GREEK SYMBOLS WILL BE RINT ED FOR ANY LETTERS ' 


state y1 


WHICH TYPE OF VARIAEBIE DO YOU WISH TO PLOT AS CIRVE NUMBER 2? 


1. STATE VARIABLE (IE., X1, X2, oftC) 
2: FEEDBACK CCNTRCI. (rf. 
32 CCNTROL INPUT (IE., U1, U2, Ere) } 
42 STATE ESTIMATE Wobsenves) {IE., XHAT1, XHAT2, FTC.) 
5. STATE BECONS TRO TION ERROR (IE., X1-XHAT1, 
X2-KHAT2, 
7 ENTER 1,2,3,4 OR 5 
? 
WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT : 
e YOU WANI TC PLOT AS THE NUMBER 2 CURVE VS TIME? 
3 
WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 
KCTE: 1. 40 CHARACTERS MAX LENGTH : 
2: GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
E) => BETA 
Fj) => PHI 
C) => THETA 
state y2 
WHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CJKVE NUMBER 3? 
1. STATE VARIABLE (IE., X1 ae cREY) 
2: FEEDBACK CONTRCL tf. tS x) 
3. CCNIRCL INPUT “) 
4. STATE ESTIMATE vdG8 gen ve Ry aoa Bre HAT XHAT2, FTC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-Xfa71 
X2-XHAT2, ETC) 
: ENTER 1,2,3,4 OB 5 
WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
‘ YOU WANT TC PLOT AS THE NUMBER 3 CURVE VS TIME? 
: 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBCLS WILL BE PRINTFD FORK ANY LSIIERS 
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BS eng eth Meta Sg ete eee Log Nie erg unig wee aie ere wee 


ENCLOSED IN eo es 


IE. (A PHA 
E 
F 
Q 


TA 
WHICH TYPE OF VARIABIE DO YOU WISH TO PLOT AS CIRVE NUMBER 4? 


>A 
> BE 
> PHI 

> THETA ) 
state y3 


1. STATS VARIABLE (IE., X1, X2,_ETC) 
2: FEEDBACK CONTRCIL (If, & = <c* 
32 CCNTRCL INPUT . U2, ETC.) 
2 STATE ESTIMATE (OBSERVER) (If., XHAT1, XHAI2, ETC.) 
5. STATE RECONSTRUCTION ERFOR (IE., X1-XHAT1, 
X¥2-XHAT2, ETC) 
= ENTER 1,2,3,4 OR 5 
d , 
WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
x YOU WANT TC PLOT AS THE NUMBER 4 CURVE VS TIME? 
2 
< WHAT IS THE CURVE LABEL FOR THIS VARIABLE? p 
NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
FE) => BETA ents 
E) => BHT > 
¢) => THETA 
state y4 
YCU MAY USE UP TO 3 HEADINGS. 
HOY MANY HEADINGS DO YOU DESIRE OW THIS GEAPH? 
Z 0, 1, 2 OR 3 , 
3 3 
WHAT IS THE DESIRED HEADING NUMBER 1? 
NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. » 
TE. (A) => ALPHA 
P} => BETA 
F) => PHI 
Q) => THETA 
closed lcorp systen 
; PHAT IS THE DESIRED HEADING NUMBER 2? ; 
NCTE: 1. 40 CHARACTERS MAX LENGTH 
2: GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
E) => BFTA 
F) => PHI 
¢) => THETA . 
example 3 
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WHAT IS THE DESIRED HEADING NUMBER 3? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 


ENCLOSED IN PARENTHESES. , 
IE. A) => ALPHA 
; => BETA 
F) => PH 
¢) => THETA 
modern control eco ry 
>> USING A FRE-ALLOCATED DATASET FOR UNIT FT17F001. 
THE FCLLCYING OPTICNS ARE AVAILABLE. : 
1. BEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
2. RKREFLOT PREVIOUSLY SAVED GRAPH DATA. 
3. EDII THE CURRENT GRAPH. 
4. PLCT REVISED GRAPH ON THE TEK618. 
5. gurt AND/OR MAKE METAFILE OF THE CUEVES. 
REVICUSLY SAVED. , 
SELECT A NUMBER EETWEEN 1 AND 5. 
5 
CO YOU PANT TO SAVE THE CURRENT GRAPE DATA IO 
EF USED LATER TO GENERATE A METAFILE? 
) 
Y OF N 
NOTE: aA METAFILE IS Be puadee FOR SMOCTH VERSATEC PLOTS. 
THERE WILL BE AN OPPORTUNITY TO GENEFATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 
y 
' 
WHAT FILE NAME DO YOU WANT THE CURVE DATA STOKED UNDER? ; 
(@ CHARACTERS MAX) 
closedlp , 
THE CORVE CATA IS BEING FILED UNDER CLOSEDIP DATA é 
END CF DISSFLA 9.0 -=- 16300 VECTORS GENERATED IN 1 PLOT FRAMES e 
FROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGC, CA. : 
1883 VIRTUAL STORAGE REFERENCES; 6 READS; O WRITES. ’ 
THE FOLLOWING OPTIONS ARE AVAILABLE: 
hs MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 
2° QUIT. 
2 ENTER 1 OR 2 
| 
WEAT FILE NAME IS THE DATA STORED UNDER? 
closedlip 
THe CURVE DATA IS BEING LOADED FROM FILE CLOSEDLE DATA 
>> USING A PRE-ALLOCATED DATASET FOR JNIT FT1T8FC001. 


TEE FOLLOWING OPTIONS ARE AVAILABLE: 
1. naAne METAFILE OF PREVIOUSLY SAVED CURVE. 


ENTER 1 0R 2 


2 
END CF DISSFLA 9.0 -- 16260 VECTORS GENERATED IN_1 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 
1886 VIRTUAL STORAGE REFERENCES; 5 READS; 0 WRITES. 
TASD 121 DETACHED 
DASD 122 DETACHED 
DASD 120 DETACHED 
aR tok ROR OR doi tok i a ok 
DC_YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
THE DISS 


PLA METAFILE THAT YOU JUST CREAT ED? 
(Y OR N 


t $128) ee 

DASD 001 LINKED R/O; E/W BY MVS; K/O BY O0085P 
Z $901) F/C - OS 

DASD 001 DETACHED 

CREATING NEW FILF: 

CREATING NEW FILE: 

FUN FILZ 6749 TO MVS COPY 001  NOHOLD 
LASD 120 DETACHED 


Oe Hee Ree oe oe a aie a ie He ae a a a a He ee a ee a eo ee a ek Ko kkk kk tok kk & 


YOUR GRAPH (S) CAN BE PICKED UP AT THE COMPUTER CENTER. 


THE GRAPH(S) WILL BE ADDRESSED TO "POP (USER ID)". 
AEC ROR IG IO Sf Ito ORO Rot Ri i a RRO RR a dk 


He Re eee ah ak ee eae hee ake ke tat oe ei oe he he ea a oe oe he a ea eo i oe ek ke a ook kK tok ok ok 


DO YOU WANT TO 


1. RUN OPISYSX AGAIN 
Ze RUN THE PLCT PROGRAM USING THE SAME MATRICES? 
3 Corr ANCTHER TYPE OF SYSTEM (OPEN/CLOSED) ) 


ENTER 1, 2 OR 3 
SR FR RRR RR RR RK BOR ORK RK RR ERK KOEI EEE K ERK ERK EEK RK 


3 


eH He A ee oe he ee ie ee ie Heke He Heo KK RR KK KK RK KK KEK KK ERK KK RK KH KKK KK KK 


HAVE A GOOD DAY!! 


SOTO RR I I IOR I SORIGR SOIT RO RR RIC IO RK HO aK 
EK; ooae 52415. 78 20:35:04 
reco off 
END RECORDING OF TERMINAL SESSION 
The graphical output generated by this example follicws 


as figure 3.2. 
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CLOSED LOOP SYSTEM 
EXAMPLE 3 
MODERN CONTROL THEORY 

TIME — SEC 
Closed-loop fime Response 


Figure 3.2 
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C. FILTER CLOSED-LOOP SIMULATION 


The following filter simulation was taken from [Ref. 7 
pp. 332 ~- 334). 

In its present cenfiguration, OPTSYSX program sequencing 
requires the input of a [CJ] matrix or design of an optimal 
regulator (if a[{S] matrix has been provided), pricr to 
initiating the oftimal estimator synthesis or user provided 
|K] matrix evaluaticn. In order to comply with Luilt-in 
program sequencing conventions, and circumvent prcgran 
difficulties which may not be specified in the particular 
system model, optimal filter synthesis may be accomplished 
by entering the identity matrix {I] in those program input 
sequences requiring the entry of an output cost (weighting) 
matrix. 

The full terminal session is recorded below, with user 
input at the left margin in lower case letters or numbers 
kelow each "2", 


ROGRAMS : 
s 
M 
T 


tIHn ty 


PTSYSX FORTRAN 
FCRTRA 


Pe yO 


NUideuNEe Ont 
abit 


HOAROM A 

Te tra OOO 
ta YSwWONMN DW rN 

HOMNIO> Hr 

Nn 

rg 


Lae 
a} 
us] 
mm 


" TC EXIT, ANY OTHER ENT 


He He te ae he ake oe ee ete ae ake hee fe hea a he oe he he KK KK RK KK RK KK 


YOU HAVE A DATA FILE NAMED ‘OPIMAT DATA' ON YOUR A 
SK SLY GENERATED BY THE OPTSYS 
Te THE F, G, H, GAMMA, A AND B 


YCU KOULD Li 


I 
i 
A 
F 
H ; E 
ESIKED DATA AT 


LD 
E 
M 
I 
T 
c 


Is 
SI 


mi. 


7 tf 


IF YOU TYPE (Y) ES. 
ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 


we ek ek ke kok ok ok kk kk KK KK KKK KK KKK KEK EEK EEK KK KEK KK KK 


PES SSE EEE EES ET SR EET SESE ST TTL SSE SEE SET ES SS FF ST TST STR SS ESS St SS SS 


DC_YCU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 


TO YOUR TERMINAL S(CREEN) OR TO A D(ISK) FILE? 
S OR D) 


PES SE SESE EES TSS SET ES ESE ELISE STEPS SPSL SESE TELE ST STR SS SES STS SS FS 


s 


he Fe he te Hee ee he eo ee ee ae ae ek KK a KK RK KKK EEK EEK KKK KKK KEKE KK KKK KE 


OUTPUT WIIL COME TO YOUR TERMINAL SCREEN. 


Me He eR oe He ee eke ke ek RK KKK KK KK KEKE KK KK KE KKH KKK KK KK KK 


LOADING CPTSYS.. 
EXECUTICN BEGINS 


OPTSYSX IS A CO 


MPLETELY INTERACTIVE OPTIMAL SYSTEMS 
EROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATIONS: 
XDOT = {F}#X + {G}#U + (GAM}* (W+W 0) 


MEASUREMENT EQUATION-- 


CCN 
ON 


Z = {H}*X + {D} #0 + V 
REGULATOR PERFORMANCE INDEX-- 
J = 1/2 * INTEGRAL (Y *{A}*Y + U *{B}*U) DT 


STATE FEEDEACK GAIN DEFINITION-- 


U = ~{C}*Xx 
DO YOU WISH TO CONTINUE? TYPE "YES" CR "NO", 
--DATA ENTRY-- 
ALTHOUGH OPTSYSX IS SPECIFICALLY DESIGNED T) READ 
ALL MATRIX LATA INTERACTIVELY, SEVERAL ALTERNATE 
METHODS ARE AVAILABLE TO USERS: 
METHCD 1--THE "F"."G" AND "GAMMA" MATRICES 
MAY BE READ PROM SEPARATE DATA FILES. 
METHOD 2--THE "F","G", AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEFINED’ WITHIN SUBROUTINE "SETUP", 
{NOTE: IN EITEEE F, THE USER SHOULD OBTAIN A COPY 
OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/R "SETIP".} 
DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO". 


DO YCU WISH TO INPUT THE "FM", "G", AND "GAMMA" 
MATRICES FrCM SUBROUTINE "SETUP" IAN THE 
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METHOD DESCRIBED ON THE PREVIOUS SCREEN? 
TYPE "YES" OR "NO", 


a GENERAL OPTSYSX GPTIONS: 
OFTION 1 -- SYSTEM puatieto WITHOUT 
CPEN-LOOP EIGENSYSTEM CALCULATIONS. 
CFTION 2 -- SYSTEM ANALYSIS WITH OPEN-LOOP 
FIGENSYSTEM CALCULATIONS. 
z CFTION 3 -- CPEN-LOOP EIGENSYSTEM FOUND 
AND PRCGHRAM TERMINATES. 
{"F"-MATRIX ENTRY FOLLOWS IMMEDIATELY. } 
CETION 4 -= MODAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OR Peouraton SYNTHESIS 
CR STEADY-STATE ANALYSIS 
el SELECT AN OPTION: 1,2,3, OR 4. 


DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
TYPE "YES" CR "NO", 


» n 
OPTSYSX LQOR/CLASSICAL OPTIONS: 
CETION 1 --— CPTIMAL FILTER AND/OR REGULATO 
SYNTHESIS WITH NO EXTERNAL "CH FOR mre 
MATRIX INPUT. 
2 CETION 2 -- CPTIMAL FILTER ae REGULATOR 
SYNTHESIS WI EXTERNAL "Cc" 
MATRIX INPUT. 
CFIION 3 -- OPTIMAL FILTER ae REGULATOR 
SYNTHESIS WITH EXTERNAL "kK" 
MATRIX INPUT. 
i CETION 4% -— CPTIMAL FILTER sae REGULA TOR 
SYNTHESIS WITH EXTERNAL "C™ AND "Kt 
MATRIX INPUT. 
7” SELECT AN OPTION: 1, 2, 3, OR 4. 
i 
> 


CC YOU WISH TO CETERMINE THE STEADY-STATE RESPCNSE 
FOR A CONSTANT DISTURBANCE? 


TYPE "YES" CR "NO", 


[TO YOU WISE TO LCETERMINE THE MODAL DISTRIBUI ION 
i] AND GAIN MATRICES? 


TYPE "YES" CR "NO", 


OPEN-LOC? TRANSFER FUNCTION OPTIONS: 
> CPTION 1 -- NO OPEN-LOOP TRANSFER FUNCTIONS COMFUTE®. 
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tom. 7 


* 


BY . 


ade.) 


~ Oe Se TO a eae ee ee De re et oe SS ee Se ee owe ee ee 


CPTION 2 -- FOLES, RESIDUES, AND ZEROS COMPUTED. 
CPTION 3 -- CNLY POLES AND ZEROS COMPUTED. 
CPTION 4 -= CNLY POLES AND RESIDUES COMPUTED. 


SELECT AN OPTION: 1, 2, 3, OR 4. 


NCISE TRANSFER FUNCTION CPTIONS: 


CFTION 1 -- NC NOISE TRANSFER FUNCTIONS 2 OMPUTED. 
CPTION 2 -- FCLES, RESIDUES, AND ZEROS COMPUTED. 
CPTION 3 -- CNLY POLES AND ZEROS COMPUTED. 

CPTION 4 -- CNLY POLES AND RESIDUES COMPUTED. 


SELECT AN OPTION: 1, 2, 3, OR 4. 


COMPENSATOR TRANSFER FUNCTION OPTIONS: 


CPTION 1 -=- NC COMP. TRANSFER FUNCTIONS 2 OMPUTED. 
CPTION 2 -- FOLES, RESIDUES, AND ZEROS COMPUTED. 
CPTION 3 -- CNLY POLES AND ZEROS COMPUTED. 
CPTION 4 -- CNLY POLES AND RESIDUES COMPUTED. 
{NOTE: A COMPENSATOR TRANSFER FUNCTION CAN BE 
CCMPUTED ONLY IF BOTH A REGULATOR 
at ILTER ARE SYNTHESIZED 


SELECT AN OPTION: 1, 2, 3, OR 4. 


WILL A FEED-FORWARD DISTRIBUTION MATRIX 
{"D" —- MATRIX} EE INPUT ? 


TYPE "YES" OR "NO". 


| 


tt Te 
mom fost hse) 


RS Sa hitdhH 
e nite 


TIFFERENT MA 
GER VALUE. 


THE DEFAULT VALUE, TYPE "0" {ZERO} 
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ob 


ve te ow ee ee De te 


DO YCU DESIRE TO SYNTHESIZE A STABLE FILTER {GCE KEGULATOR} BY 
DESTABILIZING THE ORIGINAL SYSTEM? 


{NOTE:WORKS FCK FILTER OR REGULATOR BUT NOT FCR BOTH 
IN THE SAME RUN.} 


TYPE "YES" CR "NO". 


n 
[LC YOU DESIRE TC PRINT THE ZEULER-LAGRANGE FIGENSYSTEM 
FRICR TO DECCMPOSITION {FOR CHECKING THE PROGRAM}? 
TYPE "YES" CR "NO", 
Nn 
POWER SPECTRAL DENSITY {PSD} OPTION 1 : 
OPTICN 1 -- COMPUTE THE PSD OF THE OUTPUTS AND/)R THE 
CONTROLS OF THE CONTROLLED SYSTEM FHEN FORCED BY 
PROCESS AND MEASUREMENT NOISE. (NOTE: BOTH a 
REGULATCR AND A FILTER MUST BE RESIDENT IN THE 
PROGRAM 10 USE THIS OPTION. } 
CPTICON 2 -- SAME AS CPTION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES Of EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 
OPTICN 3 -- NOT DESIRED. 
r SELECT AN OPTION: 1, 2, OR 3. 
3 
DO YOU DESIRE REGULATOR SYNTHESIS ONLY? 
TYPE "YES" CR "NOM, 
nN 
ENTER THE # OF STATES {NS} OF THE SYSTEM MAI RIX 
5 {"FN-MATRIX}. 
z 
ENTER THE # OF CONTROLS {NC} OF THE CONTROL SYSTEM MODEL 
, {"G"-MATEIX). 
4 
ENTER THE # OF MEASUREMENTS OR OBSERVATIONS {NO} OF THE 
si {"H"-MATPIX}. 
1 
ENIER THE # OF EROCESS NCISE SOURCES {NG} OF THE 
" {"GAMMA™-MATRIX}. 
1 


FIAG/FAFAMETEPR SETTINGS FOR THIS RUN ARE AS FOLLCKS: 
IOL I¢ IR ISS IM ITF1 ITF2 ITF3 IFDFH IE IDERB'G 
0 0 2 0 0 0 ) 0 ¢] 0 C 
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ISET IDSTAEB IPSD IYU INORM IREG NS NC NOB NG 


0 0 0 0 0 0 2 1 1 ; 
-{ ' 
‘ ORDER CF SYSTEM = 2 

NUMBER CF CONTRCLS = 1 

NUMBEF CF OBSERVATICNS = 1 

NUMBER CF EROCESS NCISE SOURCES = 1 


ENTER THE SYSTEM MATRIX {"F"-MATRIX} 


DIMENSION = # STATES {NS} X # STATES {NS} 
THE ELEMENT F( 1, 1)= 


2 
; 0 
J THE ELEMENT F( 1, 2)= 
| 
- THE ELEMENT F( 2, 1)= 
0 
THE ELEMENT F( 2, 2)= 
e ? 
-4.6 
THE SYSTEM MATRIX {"F"-MATRIX}... 
0.0 1:00000 
0:0 -4260000 : 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
n 
CFEN LOCP DYNAMICS MATRIX. 22. cw ccc ccc cen ecFae 
0.0 0.1000D+01 
0:0 - 0246 00D+01 
ENTER THE MEASUREMENT SCALING MATRIX [("HY-MATRIX}. ; 
le DIMEN ea = # OBSERVATIONS {NO} X # STATES {NS} 
F THE ELEMENT H( 1, 1)= 
THE ELEMENT H( 1, 2)= 
0 = 
od TEE MEASUREMENT SCALING MATRIX {"H"-MAT RIX}... | 
1.99000 0.0 ay 
4 
DO YOU FISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? “a 
TYPE "YES" CR "NO", | 
n 7 
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ere ee 


MEASUREMENT SCALING MAT RIX....22.- eee eee Hee 
0. 10COL+01 0.0 


row 
=>: 


ENTER THE OUTPUT MEASUREMENT COST MATRIX {"A"-MATRIX}. 
LIMENSION = # OESERVAT IONS {NO} X # OBSERVATIONS {NO 
THE ELEMENT A( 1, 1)= 


ae F. 
may) 


THE OUTPUT MEASUREMENT CCST MATRIX {"A"-MATR IX}... 
1.90000 


CO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


TYPE "YES" CR "NO". 


CUTEUT COST UNTER Ti ecwwious oud nds bocce eeuwes ehes 
0. 10COD+01 
ENTER THE CCNTRCL DISTRIBUTION MATRIX {"G"-MATPIX}. 
DIMENSION = # STATES {NS} X # CONTROLS {NC} 
THE ELEMENT G( 1, 1)=_ 


THE ELEMENT G( 2, 1)= 


THE CONTROL DISTRIBUTION MATRIX {"G"-MATRIX}... 
0.0 
0.78700 
DC YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


TYPE "YES" CR "NO". 


BATFIX, (uBn MAT Bm) 
ONTROLS {NC 


THE CONTROL COST MATRIX... cece ee Bowe 
1.00000 


DO YOU WISH TO CEANGE THE VALUE OF AKY MATRIX ELEMENT? 
TYPE "YES" OR "NO", 
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THE CONTROL DISTFIBUTION MATRIX..........-G.. 
9.0 
0.7870D+00 


THE CONTROL COST MATRIX ccc ccc wc cs ccces cee Bee 
0. 10COD+01 


EIGENSYSTEM OF OPTIMAL REGULATOR... ....0- 


i C-LOCE CPTIMAL REG. E-VALUES...DET(SI~F+G*C) .. ; 
-1.712C6L-01:-4.59681D+00: 
C-LOCP RIGHT EIGENVECTOR MATRIX. cc cece cca Meuee 
Pp 856588D-01 -2.125703D-01 ' 
5687503D-01 9:777458D-01 
CCNTROL EIGENVECTICR MATRIX... ccccccce ec uC ¥M.. 
-9.4S€31SD-01 3.667155D-03 
@ ' 
C-LCCP OFT. REG. LEFT E-VECTOR MATRIX..M-INV.. 
1.052798D+#00 2.292453D-01 
1:815879D-01 1:062S579D+00 
7 THE CPTIMAL FEEDBACK GAIN CONTROL MATRIX...C=BINV*GT*S... Sine 
-1.00CCD+00 -2.13495-01 
TEE CLOSED LOOP DYNAMICS MATRIX ......F-G¥*C.. 
: 0.0 1.0000C0D+00 ico 
> ~7.870000D-01 -4.768018D+00 ae 
ENTER THE PROCESS NOISE DISTRIBUTION 
MATEIX {"GAMMA"-MATRIX}. 
DIMENSICN = # STATES in} X # PROCESS NOISE SOJRCES {NG} 
. THE ELEMENT GAM( 1, 1) = 
° 0 
: THE ELEMENT GAM( 2, 1) = 
0.1 
THE PROCESS NOISE DISTRIBUTION KATRIX 
; ("GAMMA"=MATRIX}... 
] 
0.0 
6210000 
DO YCU WISH TO CHANGE THE VALUE OF AFY MATRIX ELEMENT? 
TYPE "YES" OR "NO", 
] n 
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R THE PROCESS NOISE PSD WEIGHTING MATRIX 
‘MATRIX} . 


sh 


THE PROCESS NOISE WEIGHTING MATRIX.....Q.. 
1€.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO". 


FROCESS NOISE DISTRIBUTION MATRIX.....GAMMA.. 
0.0 
0. 10C0r+00 


FCWER SPECTRAL DENSITY - PROCESS NOISE....Q.. 
9. 1000D+02 
ENTER TEE MEASUREMENT NOISE DISTRIBUTION MATRIX {"R"MATRIX}. 


DIMENSION = # OESERVATIONS {NO} X # OBSERVATIONS {NO} 
THE ELEMENT R( 1, 1)= 


0.0000001 
THE MEASUREMENT NOISE DISTRIBUTION MATRIX. ....Ree. 


0.00000 . 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 


FCWEK SPECTRAL DENSITY-MEASUREMENT NOISE..R.. 
0.10C0l-06 
ENTER THE FEEDBACK GAIN ESTIMATOR MATPIX {"K"=MATRIX}. 


DIMENSION = # STATES {NS} X # OBSERVATIONS {NO}. 
THE ELEMENT K( 1, 1)= 


2 
$5.4 
: THE ELEMENT K( 2, 1)= 
4561 


THE FEEDBACK GAIN ESTIMATOR MATRIX {"K"-MATRIX} 
S£.40000 
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4561. 00000 


DO YOU WISH TO CHANGE THE VALOE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO", ’ 


E CLOSED LOOP FILTER DYNAMICS MATRIX IS.... 


00D+01 1.0000C0D+00 
00D+03 -4.600000D+00 


EIGENSYSTEM OF OFTIMAL ESTIMATOR... ..-e-. 
C-LOCP SUBOPT. EST. E-VALUES...DET(SI-F+K*H).. 
-5.00000Dt01, 4.99984D+01: 


C-LCCP RIGHT EIGENVECTOR MATRIX... 2c 2c oo Mane c 
9.953957D-03 -1.096216D-02 
1.00000CD+00 0.0 


MEASUREMENT EIGENVECTOR MATRIX.....H (BAR) *M.. 
9.953957D-03 -1.096216D-02 


EFT E-VECTOR MATRIX..M-INV.. 


KMS VALUES OF THE ESTIMATICN ERROR... cc ccc eee 
22674042D-03 1.0868C8D-01 


DO YOU WISH TO OBTAIN A TIM 
OF THE SYSTEM YOU ARE EVA 
(Y OR N) 


NCTE: YOU MU 


ST BE LOGGED ON AT A DUAL SCREEN 
(TEK 6 18) 


TERMINAL TO UTILIZE THIS MODE. 


), G (CONTROL), H (OBSERVABLES), GAM (NCISE), 
Cost) AND B CONTROL COST) MATRICES WILL BE 
OR REENTRY TO THE MAIN OPTSYS PROGRAM. 
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Rok gb ok dodo tok ok dk ak ok ok ok 
IF YOU ARE DISSATISFIED WITH THE RESULTS 


THUS FAR ANI WOULD LIKE TO EXIT TO CMS, 
“TYPE 'y* TO EXIT- 


OTHER INPUT TO CONTI UE) 


RHR KR KKK KKH KEK EK KEKE KR KEK ERE KEKE We Me He he ae he ea eK KK kk kK 


NY 
oy 
LGADING OPTCALC. 
ZXECUTICN BEGINS 


DURING THIS SECTICN OF THE PROGRAM YOU WILL: 


- SELECT THE TYPE CF SYSTEM RESPONSE TO PLOT 
OFEN LOOP, CLOSET ae ee FILTER/REG ULATOR) 


- FOVIDE START ANTI STOP ME FOR PLOTTING CALCULATIONS 
- SELECT THE TYPE CF DRIVING eee arth aoe OR ae 
- FRCOVIDE START AND STOP TIMES FO R NG FUNCTION (S) 
- PROVIDE DRIVING FUNCTION MAGNITUDE(S). 
CLEAR THE SCKEEN TOC CONTINUE 
THE F MATRIX 
0.0 1.00000 
0:0 -4 260000 
THE G MATRIX 
0.0 : > ios 
0.78700 re 
THE C MATRIX 
- 1.00000 -0.21349 
THE BE MATRIX 
1.90000 0.0 
THE K MATRIX 
S£.40000 
4561. 00000 
THE OLLOWINS PLOTTING OFTIONS ARE AVAILABLE IF THE 
REQ 


F 
UIRED MATRICES WEFE CALCULATED IN OPTSYSX: 
1. N LOOP TIME RESPONSE 


PE 
XDCT = (F}*X + (G} *UC 

2. CLOSED LOOP TINE RESPONSE 
XDCT = [F-GeC}*X +# (S}*0C, Uv = [C}*xX 

3. OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 
XDOT = {EL #X + {6} #0C, By EX 
XHDOT = {F}*XH #{G}U 4+ {K}*{Z - Ay H} 

4. FILTER + REGULATOR CLOSED LOOP SYSTEM RESPONSE. 
XDCZ = _(FtG*C eX + (G} *0C Z = tH sk 
XHDOT = (F}*XH +(G}U + {K[*(Z - H*XH}, UD = [C} *XH 
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SELECT 1, 2, 3 OR 4. 


‘ 3 
THE (K*H) MATRIX : 
95.40000 0.0 ‘ 
45€1. 60000 20 ee 
THE CCMBINED SYSTEM F MATRIX (2*NS X 2*NS) ce 
0.0 1.00000 0.0 0.0 
0:0 -4.60000 0.0 20 
95240000 0:0 -95:240000 1.00900 
4561. 00000 0:0 -4 561. 00000 ~4.260000 
i THE AUGMENTED G MATRIX (2*NS X NC) 
0. 0 
0.78700 
0:0 
0. 7€700 
e AT WHAT TIME DO YCU WANT TO START = a 
THE TIME RESPONSE CALCULATIONS? 
F INPUT START TIME IN SECONDS. (NORMALLY 0.9) 
0 
AT WHAT TIME DO YOU WANT TO STOP ee, 
THE TIME RESPONSE CALCULATIONS? uae 
INPUT STCP TINE IN SECONDS. ey 
0 e 3 ae 
THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAV oar 
JUST SPECIFIED INTO UP TO 500 SHALL INTERVALS FOR 
TEE INTEGRATION AND PLOTTING ROUTINES. IN ORDEF 
TC SAVE COMEUTER THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED Ide "eouEe Foss IN CURVE FIDELITY. 
: HOW MANY PCINTS DO YOU WANT TO CALCULATE? 
500 
DOES THE SYSTEM UTI1IZE A DRIVING FUNCTION (CONTROL INPUT)? 
: (Y)ES OR (N)O 
: THO TYPES . cE FUNCTIONS CAN BE USED AS DRIVERS. a 
le 1. STEP INPUT 
ra BARE INP DT 48 
ENTER YOUR SELECTION, 1 GR 2. FOR DRIVING FUNCTION NUMBER 1 con 
1 


AT WHAT TIME DO YOU DESIRE INPUT NUMBER 1 T9 START? ; 
INPUT THE START TIME 1N SECONDS. 
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- THE STATES, 
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WITH THE RESULTS 


THUS FAR ANT WOULD LIKE TO EXiTt TO CAS, 


TO EXIT- 
THE PROGRAM PLOTS 


yt 
(ANY OTHER INPUT TO CONTINUE) 
FO goo kkk algo gg ak fk dak dk doo kk kk dk ke 
execution will soon follow ... 


-TYPE 
F 
R 
EXECUTICN BEGINS... 


YOU ARE DISSATISFIED 
«+. Your Fortran program is now being loaded ... 


THIS PORTICN O 


7 
iF 


NMOS DONO ST NMNM = PNOWE[FDAVHNOOSTNAAMM SWMNOOlPDONVAGD 
DQDQNQ0 OO CC CF FF TBE ME NAN NANNANNNANNNANCINGS 
ee eeeeeeertsnstst err? eeeseeeertrteest eet te ee eeeerseseeeeet 8 @ 


Se She eae a ok oe ee ote ate ake feo eo ae aa ee ae a ote 2 ake fe ic ae a a a ak ie ake ake ak ack oe oe ake fe i ake ake ake ic ok oe 


S ALSO AVAILABLE TO KEVIEW ANY 
wea PREVIOUSLY SAVED AS DATA 


CLEAR THE SCREEN TO CONTINUE. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. PLCI THE DATA YOU JUST CALCULA 
2- PLCT A CURVE THAT YOU PREVIOUS 


ENTER 1 OR 2 


TED. 
LY SAVED. 


YOU MAY PLCT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DO YOU WISH TO PLOT? 


n 


a 


WHICH TYPE OF VARIAFPLE DO YOU WISH TO PLOT AS CURVE NUMBER 1? 


VARIABLE (IE., X1, X2, ETC) 
CK CONTRCL (IE. = =¢# 
L INPUT U2 


TIM (OBSER VER) XHAT2, ETC.) 
CONSTRUCTION ERKOR 
es 


TC) 
TER 1,2,3,4 OR 5 


WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TC PLOT AS THE NUMBER 1 CURVE VS TIME? 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


1. 40 CHARACTERS MAX LENGTH 
2. GREEK SY BE PRINTED FOR ANY LETTERS 
ey D ESES. 


_ 


angular position - 

LE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 
pee cele 

4, 02, 
R 


tt 


WHICH 


Hi 


ETC.) 
vip) (IE., XHAT1, XHAI2, ETC.) 
ERR (IE., Xt-xfaT1, 


<M Mh 
NOHMAHIOMH K 
team nit  ont 
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if 


( 


Bm). 


ade.) 


ENTER 1,2,3,4 OR 5 


WHAT IS THE SUBSCRIFT OF THE STATE ESTIMATE THAT 
YOU WANT TC PLOT AS THE NUMBER 2 CURVE VS TIME? 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
eT IN Eee ieree . 


=> ALP 
E) => BETA 
F) => PHI 
wat ¢) => THETA 
angular frosition estimate - (x)e1 


YOU MAY USE UP TO 3 HEADINGS. 
HOW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 


0, 1, 2 OR 3 


3 
WHAT IS THE DESIRED HEADING NUMBER 1? 
NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. {&) => ALPHA 
E) => BETA 
F) => PHI 
Q) => THETA 
filter only closed lcop 


WHAT IS THE DESIRED HEADING NUMBER 2? 


NOTE: 1. ae CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED ete 


IE. > 
ay 53 
¢) => 


+ 


HI 
HETA 


B 
P 
T 


example 4.1 


WHAT IS THE DESIFED HEADING NUMBER 3? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. . 
IE. (A) => ALPHA 
E) => BETA 
F) => PHI 
Q) => THETA 
r_optimal contrcl systems 
ING A PFE-ALLOCATED DATASET FOR UNIT FT17F0) 1. 
E FCLLOWING OPTICNS ARE AVAILABLE. 
EEGIN NEW GRAPH CF OTHER CONTROLS STATES, OR ESTIMATES. 
REFLCT PREVIOUSLY SAVED GRAPH DATA. 
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Ss ww we es 
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OP ESTIMATES. 


NTROLS, STATES, 
PA DATA. 
ON THE TEK618. 


ETAFILE OF THE CURVES. 


F OTHER CO 


Om ° 


TRE FOLLOWING OFTIONS ARE AVAILABLE. 


SELECT A NUMBER BETWEEN 1 AND 5. 


? 
3 


THE GRAPH ELIT MENU 


CURVE ON THE CURRENT PLOT. 


OQ 

e 

Ns) 
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w 
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isa) e 
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AMHARESHAOR 


He « HHS 
Ba HMOWNO 
DPI WTO 
POmMHA MEA, 
De ro) 
1 YO Cubs way = 
VOUVUH 


TwOQrmnaeaeg 
VARNUULVVLNUHN 


eeoese#ee?e#ereer#e Ff 
[AMANO AHO 
<= 


ELECT A NUMBER RETWEEN 1 AND 10. 


S) 
2 
5 


CN WHICH CURVE CC YOU WANT TO CHANGE THE Y-SCALE? 


OR 4 


3, 


ENTER CURVE NUMBER- 1, 2, 


Y-MIN VALUE AT THE ORISIN? 


WHAT IS THE NEW 


WHAT IS THE NEW Y-MAX VALUE? 


ele 


' 
a 


THE GHKAPH EDIT MENU 


CURVE ON THE CURRENT PLCT. 
.5 X 6.0) 


A 
URRENT PLOT. 


OR ADD 


Ge MMORNO 
gm > hmm aAmou 
PUmImMRiEIQeA 

Dg © 
HWO Mabe 
Or OUUOOH 
2ARHSZaaZASH 
SH HeQetetetiecH 
MNVAATMAUMNA 
YARINUVUYNUVYA 


TAM FIOM OHNO 
= 


10. 


SELECT A NUMBER PETWEEN 1 AND 


10 


THE FCLLOWING OPTIONS ARE AVAILABLE. 


OR ESTIMATES. 


STATES, 


LE OF THE CURVES. 


HE TEK618. 
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FREVIOUSLY SAVED. 
SELECT A NUMBER PETWEEN 1 AND 5. 
4 
E FOLLOWING OPTIONS ARE AVAILABLE. 


N NEW GRAPH OF OTHER CONTROLS, STATES, Ox ESTIMATES. 
C Pe ceee ena na D GRAPH DAT. 


ON THE TEK618. 
ETAFILE OF THE CURVES. 


7 
_., JUST B 
filteron 


YOUR ANSWER MUST BE "¥" OR "NM, 


DO YOU WANT TO SAVE THE CURRENT G 
BE USED LATER TO GENERATE A METAF 


Y OR N eee 


WHAT FILE NAME CC YOU WANT THE CURVE DATA STOKED UNDER? ere 
(8 CHARACTERS MAX) 


filteron 
THE CURVE CATA IS BEING FILED UNDER FILTERCN DATA 
END CF DISSELA 9.0 -- 43644 VECTORS GENERATED IN 3 PLOT FEAMES 
ERCPRIEIARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 
8493 VIRIUAL STORAGE REFEPENCES; 9 READS; O WRITES. 


THE FOLLOWING OPTIONS AKE AVAILABLE: 


1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 
2. QUIT. 


ENTER 1 OR 2 


: WHAT FILE NAME IS THE DATA STORED UNDER? 
filteron 


THE CURVE DATA IS 
>> JSING A PRE-ALL 


ip Pie oe a ey AS Beg tye Pe Bae OA re 
BSE gh. Pe irish yaa tees em ee ey OL Sh ee pee 
re eg et EE eS en rg 2 


FE FOLLOWING OPTIONS AEE AVAILABLE: 
MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 


§ ENTER 1 OR 2 , 


14919 VECTORS GENER 
RE PRODUCT OF ISSCO, 
EFERENCES; 9 READS; 


PLOT FRAMES 
’ 


PS SPS SESE ET ESTES STS SESE SESS EES S ESS SETS SE TTS TTS TS TR SST SS STS tS SS 


DO_YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
TEE DISSPLA ar On N) JUST CREATED? 


eR Me eK ee oe ee ee oe ae fe i Re a Re Ke ee Re ae Ra a ea Ro a ea oe ie 


4 120 ¢g 
DASD 001 LINKED R/O; R/W BY MVS 
Z $901) RC - 0S 
pe DASD 001 DETACHED 
CREATING NEW FILE: ’ 
CREATING NEW FILE: 
PUN FILE 6910 TO MVS COPY 001 NOHOLD 
LASD 120 DETACHED 


Me a ae fe ke Re fe 2 fe ai ae te ee ae ke ae oi he he ae ah ae oe ofeach ea eo ete eo oo te oa ok ok 


YCUR GRAPE (S) CAN BE PICKED UP AT THE COMPUTER CENTER. 


THE GRAPH(S) WILL BE ADDRESSED TO “POP (USER ID)". 
3 igok kk fo dk ia dak idk ioc doi iki ig lok kai to doko tok dato a ok 


Ke He He He Ree oe ee ee ee ae a KH a ok toe ee ke 


DO YOU KANT TO 


1. RUN OPTSYSX AGAIN 
2. RUN THE FLOT PROGRAM USING THE SAME MATRICES? 
(TO_ PLOT ANOTHER TYPE OF SYSTEM (OPEN/CLOSED) ) 
. 3. QUIT 
ENTER 1, 2 OR 3 

S FC RIG RRR II IR RR IRR SOR RO ok do ka 
2 3 
. 3 IIR fog tok ko Sito fk OR tok Rok tk kok ok kk 


HAVE A GOOD DAY!! 


fete SS SLES Pe TTS SSS SSS ETE SLSR SESS LS SS SSS ES SS STS SS TSS SLES SESS FS Ss 
Rs £=13.37/21.24 21325309 

record off 

END RECCFDING OF TERMINAL SESSION 
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The graphical output generated by this example follcws 


as figure 3.3. 


FILTER ONLY CLOSED LOOP 


EXAMPLE 4.1 
LINEAR OPTIMAL CONTROL SYSTEMS 


ost'o 


ogto 


Scro 


get'o 


0010 $200 060°C $z0'0 000°0 S20'0~ 050°09240'0- 


1? -— NOILISOd UVYTNONYV 


Oe ae OE ee Ee ee ee ee gars Se ee ten ee, 
Oot'o 920 0 0¢0°9 9200 oac a $20'0- 090°04940'0- 
Td? — ALVWILST NOILISOd UVIADNV 


0.16 0.20 0.26 0.30 
TIME — SEC 


-10 


0.06 


0.00 


Figure 3.3 


Pilter Closed-loop Time Response 


D. FILTER AND REGULATOR CLOSED LOOP TIME RESPONSE 


The following filter and regulator example was taken 
from [Ref. 7 pp. 382 - 384]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case letters or numbers 


below each "2", 


record on 
BEGIN RECORDING OF TERMINAL SESSION 
Fs T=0.01/0.02 09:08:30 
optsys 
THE OFTSYS EXEC CONTROLS A TRIO OF PROGRAMS: 
1. OPTSYSX FCRTRAN SYSTEM Anais iS). 
2. OPTCALC FCRTRAN CALCULATE TIME eae 
3. OPTPLOT FCRIRAN DISSPLA PLOTTING ROUTINE) 
EACH PRCGRAM PASSES INFORMATION TO THE NEXT PROGEAM 2 
THROUGH A DATA FILE WRITTEN TO THE USERS DISK. IN TAIS 
CASE, THESE FILES ARE “"OPTMAT DATA" AND “OPTPLOT DATA". 
THE SIZE OF THESE FILES VARY WITH THE SYSTEM ORDZR, AND 
CAN USE ABOUT 20% OF THE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 
- TYPE "E" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 


RA He Ke Be Re Re ee eae ee Heke ee eK eK ea a a ek ete eK 


YOU HAVE A DATA FILE NAMED ‘OPTMAT DATA* ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
FROGKAM AND CCNTAINS THE F, G, H, GAMMA, A AND 3B 
MATRICES FROM THAT RUN. 


IF YCU WOULD LIKE TO USE THESE SAME 

THIS RUN, THE  CPTSYS PROGRAM WILL 

LESIRED DATA AT THE APPROPRIATE TIME, 
IF YOU TYPE (Y) ES. 

ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 


PES SES PSS ESS TSS S STS SET ET TTS STEP ET SESE STS TTS TTS TT SSS ES SS EES tS SS 
bi 


HATRICES 
READ I 


oe He eee he he fe ie he eo ate ee ie a ee a ee he oe a ee a ae ea a a ee he ee ek ake ke 


DO YOU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 
TO YOUR "ERMINAL 34 °RE*5} OR TO A D({ISK) FILE? 


PES S SESE STE ST TSS ECE ESTEE TIS ESSE SLES ESS TES TSS LST ES ESAT SSS tS tS SS 


s 


HK KR HK AK HK KK RR KK RR RK KK RR RR KKK KK eK KE ER RK KK 
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Ete as. See SS ee a, Eee top at kh ee eae RSE 8) Rng eh gee NG a a Pat ae ihe at) ERR eae asd Puig pe Ram RY Sie, Peg) aay es | leh tags eee Ooh eh 
AP. BAe Bretagne ie tet Gee Dating race hg ORG Ti, Sete pay, Bi een er, Wee Oem wig hO owe ta. 8 ete wre er, We Cree Pan aig tome aay chine Pek eee ane ieee meg), Suey: BERS, BN eee. ao 


OUTPUT WIIL COME TO YOUR TERMINAL SCREEN. 
go kk took ok kt ok ok ok ok ok kk 


LOADING CPISYS 


EXECUTION BEGINS... 
OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CCNTROL 
EFROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES CF ey eEne CONTROL EQUATIONS: 
XDOT = [{F}*X + {G}*U + {GAM} *(W+W9) 
MEASUREMENT EQUATION-~ ' 
2 = {H}*X + {D}*0 + V 
REGULATOR PERFORMANCE INDEX-- 
J = 1/2 * INTEGRAL (Y *{A}*Y + U *{B} *0) DT 
STATE FEEDBACK GAIN DEFINITION-- 1 
U = -{C}*X 
DO YOU WISH TO CONTINUE? ‘TYPE "YES" CR "NON, 
Y --DATA ENTRY-~ 
ALTHOUGH OPTSYSX IS SPECIFICALLY DESIGNED TO READ 
ALL MATRIX DATA INTERACTIVELY, SEVERAL ALTERNATE 
METHODS ARE AVAILABLE TO USERS: 
METEOD 1--THE "E","G") AND "GAMMA" MATRICES 
MAY BE READ FROM SEPARATE DATA FILES. 
METHOD 2<-THE "F","G", AND "GAMMA" MATRICES MAY BE seman 
EXPLICITLY DEF{NED’ RITHIN SUBROUTINE "SETUP". 
; (NOTE: IN ETTEER CASE, THE USER SHOULD OBTAIN A COPY 
OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/R "SETUP".} 
DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO", 
y sas 
pO YOU FISH TO INPUT THE "F", "G", AND "GAMMA" 
MATRICES FRCM SUBROUTINE "SETUP" IAW THE 45 
METHOD DESCRIBED CN THE PREVIOUS SCREEN? oo 
TYPE "YES" OR "NO", 4 
nN 
GENERAL OPTSYSX OPTIONS: 
OEFTION 1 -- SYSTEM ANALYSIS WITHOUT ee 
OFEN-LOOP EIGENSYSTz s éaLcuLar Ions. ae 
CETIGN 2 —- SYSTEM ANALYSIS #ITH OPEN-LOOP y 
EIGENSYSTEM CAL ULATIONS. ; 
CETION 3 — OPEN-LOOP EIGENSYSTEM FOUND on 
AND PROGRAM TERMINATES. EY 
{"F"=MATRIX ENTRY FOLLOWS IMMEDIATELY. } a 
OFTICN 4 — MCIAL DISTRIBUTION MATRICES COMPUTED . 
WITHOUT FILTER OB REGULATOR SYNTHESIS ay 
OR STEADY-STATE ANALYSIS. zc at 
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eng ee oneg. thee Mg a RIN. oes Senet age Ee eR ye tn Me eRe ath 
eh sel ac) Lie petien. Stet eer Oy peg eh mts ee Gar ie Me ote Ae eB iy aps a mt ete a SD Te ER Cae rn en alte te eh ee oy 
SR aS PREC ad Rc AR a al ek Ps Rr ea te 


Ce ee eee ee 


SELECT AN OPTION: 


1,2,3, OR 4. 


DO YCU DESIRE RMS VALUES OF STATE AND CONTRIL? 


TYPE 


CETION 


OFTION 


CFTION 


OFTION 


SELECT 


1 


“YE Ss" OR "NO . 


CPTSYSX LOR/CLASSICAL OPTIONS : 


-- OPTIMAL FILTER 


AND/OR_ REGULATOR 


SYNTHESIS WITH NO SXTERNAL "C* OR "KN 
MATRIX INPUT. 

2 -- OPTIMAL FILTER AND/OR REGULATOR 
SYNTHESIS WITH EXTERNAL "C" 

MATRIX INPUT. 

3 -- OPTIMAL FILTER AND/OR REGULATOR 
SYNTH FSIS WITH EXTERNAL "K" 

MATRIX INPUT. 

4 -- OPTIMAL FILTER AND/OR REGULATOR 
SYNTHESIS WITH EXTERNAL "C" AND "KN 
MATRIX INPUT. 

AN CETION: 1, 2, 3, OR 4. 


[CG YOU WISH TO LETERMINE THE STEADY-STATE RESPONSE 
FOR A CONSTANT DISTURBANCE? 


TYPE 


AND GAIN MATRICES? 


TYPE 


wyEsSn 


CR 


"YES" CR "NO", 
DO YOU WISH TO LETERMINE THE 


tno ", 


MODAL DISTRIBUTION 


OPEN-LCCP TRANSFER FOUNCTICN OPTIONS: 


OEFTICN 1 
CETION 
CETION 
CETION 


2 
3 
4 


NO 


SELECT AN OFTION: 


OPEN-LOOP TRANSFER FUNCTIONS COMPUTED. 
POLES, KESIDUES, 
ONLY POLES AND ZEROS COMPUTED. 

CNLY POLES AND RESIDUES COMPUTED. 
14. 25°33 


AND ZEROS COMPUTED. 


OR 4. 


NOISE TRANSFER FONCTION OPTICNS;: 
NC NOISE TRANSFER FUNCTIONS COMPUTED. 


OFTICN 1 
CETION 2 
OFTICN 3 
CETICN 4 


PCLES, 


RESIDUES, 


AND ZEROS COMPUTED. 


ONIY POLES AND ZEROS COMPUTED. 


ONIY POLES AND RESIDUES 


COMPUTED. 
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a ee 


SO SR Ee eee elie Sot i a ee al 
. sae et ae ee 


a eee ae 


SELECT AN OPTICN: 1, 2, 3, OR 4. 


COMFENSATCR TRANSFER FUNCTION OPTIONS: 


CFTICN 1 -—- NO COMP. TRANSFER PUNCTIONS COMPUTED. 
CEFTICN 2 -—- PCIES, RESIDUES, AND ZZROS COMPUTED. 
OFTICN 3 -- ONLY POLES AND ZEROS COMPUTED. 

CETICN 4 -- ONIY POLES AND RESIDUES COMPUTED. 


{NOTE: A COMPENSATOR TRANSFER FUNCTION CAN EE 
CCMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
ANL/OR INPUT. } 


SELECT AN OFTION: 1, 2, 3, OR 4. 


WILL A FEED-FORWARD Bae eet TON MATRIX 
{"D" - MATRIX} EEF INPUT 
TYPE "YES" OR "NO", 
n 
THIS OFTICN DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKOV PARAMETER IS ZERO-THE MARKOV PARAMETER INDICATES 
peer en OF THE NOUMERATOR POLYNOMIAL OF EACH TRANSFER 
ALL "N" ZEROS OF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0. 
LESS THAN 10.0**{-IEF} IS CONSIDERED ZERO. 
THE DEFAULT VALUE OF THIS PARAMETER {IE} IS 6. 
IN CTHEE WORDS, IE = 1.0E-6. 
IF YOU DESIRE A DIFFERENT MARKOV CR™TERIA, 
TYPE THE INTEGER VALUE. 


IF YCU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 


DO YCU DESIRE TO SYNTHESIZE A STABLE FILTER {Oh REGULATOR} 
DESTABILIZING THE ORIGINAL SYSTEM? 


{NCTE:WORKS FCR FILTER OR REGULATOR BUT NOT FOR BOTH 
IN THE SAME RUN.} 


TYPE "YES" CR "NO". 


DO YOU DESIRE TC PRINT THF EULER-LAGRANGE EIGE 
FRICR TO DECOMPOSITION {FOR CHECKING THE PR 


TYPE "YES" CR "NO", 


POWER SPECTRAL DENSITY {PSD} OPTION 1 : 
OPTICN 1 -- COMPUTE THE PSD OF THE OUTPUTS AND/O& THE 
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ee ee ee ae ee ay 


BY 


{ 

{ 

4 

{ 
a 
; 

4 

‘ 

* 


PR a ee ne, Oe Ge ee 


Ce er Se Se Si RO a Re OR Ee DR a 


ae ean ea 


CONTROLS OF THE CONTROLLED SYSTE4 WHEN FOR 
PROCESS AND MEASUREMENT NOISE. {WOTE: BOTH A 
REGULATOR AND A FIITER MUST BE RESIDENT IN THE 
PROGRAM IC USE THIS CPTION. ‘ 
OPTION 2 -- SAME AS CETION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES CF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. - 
CPTICN 3 -- NOT DESIRED. o 
: SELECT AN CFTION: 1, 2, OR 3. 
3 
TO YOU DESIRE REGULATOR SYNTHESIS ONLY? 
TYPE "YES" CR "NO". 
n 
THE wre, nen, tye e "GAM", wan AND “ EW y ATR ICES 


FROM YOUR PREVIOUS OPTSYS RUN WERE SAVED. 


THE FOLLOWING See eS Bet AVAILABLE: 
1. USE ALL OF THE SAME MATRICES ASAIN. 


ENTER 1, 2, OR 3. 


NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 

AND YOU WIIL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


FLAG/PARAMETER SEITINGS FOR THIS RUN ARE AS FOLLOWS: 
IOL I¢ IR ISS IM ITF1 ITF2 ITF3 IFDFW IE IDEBUG 


0 0 3 0 ¢) 0 0 0 0 0 0 
ISET IDSTAE IPSD IYO INORM IKEG NS NC NOB NG 

0 0 0 0 0 0 2 1 1 1 
ORDER OF SYSTEM = 2 
NUMEER CF CONTROLS = 1 
NUMBER CF CBSERVATICNS = 1 
NUOMEER CF FROCESS NOISE SOURCES = 1 

THE SYSTEM MATRIX {"F"-MATRIX}... 
8:6 4720008 


rO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
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TYPE “YES" OR "NO", 


OPEN LOOP DYNAMICS MATRIX. cc cecce cece cece eFes 


Q.0 0.1000D+01 
0.0 -0.46005+01 


THE MEASUREMENT SCALING MATRIX {"H"-MAT RIX}... 
1.00000 0.0 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE “YES" CR "NO". 


MEASUREMENT SCALING MATRIX. ccccecc cee cece ele. 
0. 10COr+01 0.0 


THE CONTROL DISTRIBUTION MATRIX {"G"-MATRIX}... 


0.0 
0.78700 


DC YOU WISH TO CEANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO". 


ENTER THE FEEDBACK GAIN CONTROL MATRIX {"C"-MATRIX}. 


DIMENSION = # CONTROLS {NC} X # STATES {NS}. 
THE ELEMENT C{ 1, 1)= 


? 
-254 


eT 
= THE ELEMENT C( 1, 2)= 
-19.57 


THE FEEDBACK GAIN CONTROL MATRIX {"C"-MATRIX} 
-254.1000C ~ 19.57000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 


THE CONTROL DISTRIBUTION MATRIX... ...e2 400 oGee 


0.9 
0.7870D+00 


THE CPTIMAL FEEDBACK GAIN CONTROL MATRIX...C=BINV*GT*S... 
-2.5410D+02 -1.9570L5+01 


THE CLOSED LOOP D 


0.90 1.0000 
-~1.999767D+02 -2.0001 


YNAMICS MATRIX ..-..2. F-G*C.. 
00D+0 
59D+0 

C-LOCP SUBOPT. REG. E-VALUES...DET(SI-FtG*C).. 
-1.000C8D+01, 9.99804L+00: 


C-LOCP PIGHT EIGENVECTOR MATRIX... .ceee ewe Mewes 
-5. 8QD- ~4 .999602D-02 
1: 398838b+86 ~ $26 
CCNIFCL EIGENVECTCR MATRIX... ccccc cece e eo eC ¥Me. 
-6.862510D+00 1.270399D+01 


C-ICCP SUBOPT-REG. LEFT E-VECTOR MATRIX..M-INV 


9.0 1.009000D+00 
-2.00015SD+01 -1.000276D+00 


THE PROCESS NOISE DISTRIBUTION MATRIX 
{"GAMMA"—MATRIX}. oe. 

0.0 

0.10000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 


ENTER THE PROCESS NOISE PSD WEIGHTING MATRIX 
{"C'MA TRI X} . 


DIMENSION NOISE SOURCES gy? Xx 


NOISE SOURCES { 
THE ELEMENT = 


THE PROCESS NOISE WEIGHTING MATRIX.....Q.. 
10. 00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPES "YES" CR "NO", 


FROCESS NOISE DISTRIBUTION MATERIX.....GAMMA.. 


Q.0 
0. 10C0D+00 


POWER SPECTRAL DENSITY - PROCESS NOISE....Q.. 
0.10005+02 
ENTER THE MEASUREMENT NOISE DISTRIBUTION MATRIX {"R"MATRIX}. 


DIMENSION = # OESERVATIONS {NO} X # OBSERVATIONS {NO} 
THE ELEMENT K( 1, 1)= 


z 
0.9000001 


THE MEASUREMENT NOISE DISTRIBUTION MATRIX... Roe 
C.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


TYPE "YES" CR "NO". 


POWER SPECTRAL DENSITY-MEASUREMENT NOISE..R.. 
0. 10COLD~06 
ENTER THE FEEDBACK GAIN ESTIMATOR MATRIX ["K"-MATRIX}. 


DIMENSION = # STATES {NS} X # OBSERVATIONS {NO}. 
THE ELEMENT K( 1, 1)= 


rs 
95.4 


if THE ELEMENT K( 2, 1)= 
4561 
THE FEEDBACK GAIN ESTIMATOR MATRIX {"K"-MATRIX} 


9£.40000 
4561.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 


FILTER STEADY STATE GAINS. ccc cc cece ccc eco Kone 
9.-540000D+01 
4.56 10CO0D+ 03 


CLOSED LCOP FILTER DYNAMICS MATRIX IS.... 


D+01 1.000000D+00 
D+03 -4.600000D+00 


EIGENSYSTEM OF OFTIMAL ESTIMATOR... ...e0. 


C-LOCP SUBOPT. EST. E-VALUES...DET(SI-F*K*H).. 
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—5.00000L5+01, 4.99984D+01: 
C-LCOF RIGHT EIGENVECTOR MATRIX. cece ese ceMecee 


7D-03 -1.096216D-02 

CpD+00 0.0 
MEASUREMENT EIGENVECTOR MATRIX.....H (EAR) *M.. 
9.953957D-03 -1.096216D-02 


C-1ICCP SUBOPT. FIITER LEFI E-VECTOR MATRIX..M-INV.. 
0.0 1.0 
-9.122292D+C1 9.0 


COVARIANCE 


RMS VALUES OF THE ESTIMATION ERROR. ccc ceeceee 
2-67404ZD-03 1.0868C8D-01 
DO YOU WISH TO OBTAIN A TIME RESPONSE 
CF THE SYSTEM YOU ARE EVALUATING? 
{Y OR WN) 


OU E LOGGED ON AT A DUA 
EK ERMINAL TO UTILIZE T 
YS (SgBTRO 1 OBSER 
TE AND ute L COs 
VE ENTRY TC THE MAIN O 


NCTE: 


Y 
(T 


THE F & 
A (OU 
SA 


Y 


He ee oe Re oe he eRe oh oe oe ae Se hee ete He he He oh fe he eke Ee fe he ee he eR he ae ah eke a a a ea ak a oe eK oe tek a Ke 


IF YCU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WCULD LIKE TO EXIT TO CMS, 


“TYPE ‘¥* TO EXIT- 


ANY OTHER INPUT TO CONTINUE 
ete TTTTCCSt Te TT Stet tet Tet eter et TS: 


LOADING CPTCALC.... 
EXcCUTION BEGINS... 


ine deaneet ae sancanseese 
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rary 
ey? 


CLEAR THE SCREEN TO CONTINUE 


THE F MATRIX 
0.9 1.00000 
0:0 -4260000 
THE G MATRIX 
0.0 
0.78700 
THE C MATRIX 
-254.10000 -19.57000 
THE H MATRIX 
1.00000 0.0 
THE K MATRIX 
95.40000 
4561. 00000 
TRE FOLLOWING PLOTTING CPTIONS ARE AVAILABLE IF THE 
REQUIRED MATRICES WERE CALCULATED IN CPTSYSX: 
1.  OFEN LOOP TIME RESPONSE 
XDOT = {F}*X + {G}*UC 
2. CLOSED LOOP TIME RESPONSE 
XDOT = [F-G¥C}*X + {G}¥*UC, U-= {C} *X 
3. OEFTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 
XDOT = {F}*X + {G}*UC, =" {H} *X 
XHDCT = {F}*XB +{G}U + {K}*{Z - B*XH} 
4. FILTER + REGULATOR CLOSED LOOP SYSTEM RESPONSE. 
XDOT =_{FeGeC)#X + fc} * - (B}s 
RHDCT = {F}*#XH +{G}U + {Kf*{Z - H*XH}, U = [C}*XH 
A SELECT 1, 2, 3 OR 4. 
4 
THE (G#C) MATRIX 
0.0 0.0 
-1992€7670 -15.40159 
THE (K*H) MATRIX 
$8.40000 0.0 
4561. 60000 0.0 
THE COMBINED SYSTEM F MATRIX (2*NS X 2*NS) 
0.0 1.00000 0.0 0.0 
020 -4260000 -199.97670 -15.40159 
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ae i ee eae Lae gO) 
ast oot ee a ew 


en eee 


$5.40000 0.0 -95.40000 1.00000 
4561. 00000 0.0 -4760.97670  -29.00159 
r| THE AUGMENTED G MATRIX (2*NS X NC) - 
0.9 
0.78700 
0.0 
0.78700 
AT WHAT TIME DO YOU WANT TO START 
THE TIME RESPONSE CALCULATIONS? 
INPUT START TIME IN SECONDS. (NORMALLY 0.0) 
0 
‘ AT WHAT TIME DO YOU WANT TO STOP 
TEE TIME RESPONSE CALCULATIONS? 
: INPUT STCOE TIME IN SECONDS. 
0.6 
THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAVE : 
JUST SPECIFIED IKTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
TO SAVE COMPUTER TIME, THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 
7 HCW MANY POINTS DO YOU WANT TO CALCULATE? 
500 
DOES THE SYSTEM UTILIZE A DRIVING FUNCTION (CONTROL INPUT)? 
(Y) ES OR (N)O 
n 
t,] DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 
(Y)ES OR (N)0? 
nN 
: WHAT IS THE INITIAL CONDITION FOR X( 1) ? 
Ou 
z WHAT IS THE INITIAL CONDITION FOR XHAT( 1) ? 
) 
; é WHAT IS THE INITIAL CONDITION FOR X( 2) ? 
le 0 
3 WHAT IS THF INITIAL CONDITION FOR XHAT( 2) ? 
0 


THIS IS YOUR LA TO 
MAXE CHANGES IN THE FOLLOWING A&EAS. 
1. SELECT ANOTHER TYPE OF SYSTEM TO PLCT 
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ag ae ae NS Se PN eo ae GSES Faves hl ee Seat ee onret cat ete eaee net ak ee RT ORAS La ah eR ar RRA SS ao RRS a anaes | Ras Sa Aerie as 
tes ree me rang eer fea cen Te Sieg ig. Mee OL arte eras Mim Beg Sed Mas Wg Pe eee Me agin Mites ee re wT antl | CE STORE RR RE te ee 2 San ee Sa Ste 
Baa ee, iy, See we Tee tT eS Pie Beg ey ie at er ee et ae ae TS, Soa Re he ee Ce he elt ele ae eee eto oe ual Sg Sag te 


— sa = 


FILTER OR FILTER/R EGULATOR) 


(OPEN, CLOSED, 
AND STOP TIMES 


START 


2. 


DRIVING FUNCTIONS 
INITIAL CCNDITIONS 


CONTINUE 


35 


4. 


5. 
SELECT A NUMBER BETWEEN 1 AND 5. 


CATICN OF PROPER PROGRAM OPERATION. 
STATES AND STATE ESTIMATES CAN BE PLCTTED. 


NG INFORMATION IS PROVIDED INLY 


Hi 


ALL CCNTROLS, 


X (3) 


X (2) 


X (1) 


u(1) 


TIME 


ODF NNN GEMM TEE HK NNANANANAANNAMMMMMMMAMMMMMMNMM 
DODOODO DOCS VOSCOSOCOCDDVDOD SCO CDCVOVODVDOSCVOVOOVOOVOOVOOOO0OO 
t(terertrb’evrebtebtetretr;etretteteterbrtredaebprrtetrtebretet 
ARAARAARAARAARARAARARAARARARAARAARAARRAARARARAARAARAARARgQAAARAAARARAA 
DFM NAAD NM OM SH NEW DW DORMNAPENHNDO FNAVOEN NOOO Ae TOMAHAO 
DTNOMDMNONT HNODOMEE & OME KM COT ONVDDHNOSFORMSCAMDOM CMON 
NOAM OMNODE POST THA CONNNVOANATWOANANTKMNT ME ODNWODDe OOMNWO 
ONNARMONK OF TOW FAH AMMNNAANNCUIK DOP NEOPMPM ZTANHOP DE HADOHANIM 
DFHNOWNNAE KADNNFTDNOMN OM OME M QNOMNIDANNINE NONE SOCWUPKOANIHMN 
TOMTDOR. NNN NA NODDONANNKDMNNNEKNMNDHMNMNIOMNMNOK OH KNCAMAN ss 
DOPOD IMC FINE CF HME FHKE N OF ON PNAC OO TEKH NNWOE EE OOM 
eeeseeveeteseeeeoese*teeeeee#rsteeeeteeereeeeetrtereeest © Ce &« Fe Ce © eo @ 
ololelelelelelelelelelelelelelelelelelelolelelelolelelelelel ol Floleltloltlalelelelelelalelelela) 
preter er ee ever eherovaerevaeraed 
OOK KRM COO OOO CEPR OOOO CEE PEER ER EAAAMANAING 
HelelelelelelelelelelelelelelSlelelolelelelelelelaelelelelelaleleleleleleleleilelelolelelelele) 
+eeteeeeeeeee ei Pir itr +eeeetirrrr+rrirtrt)etortrrtrtrt)errertere 
AAARRAARARAARAARAARAARARAARNAARANRAARARAANARAAQAAANAANAARAARAAANAAAG 
MONTNMAMNADNFOONDODODKRMANFOWODNNAS F FODOANNAMANMNAM OM +> 
TMK DMNHDNHNOMOM FADDOM OM K- ON DODMNHOWNKHDMNMNON™ SHOMEK MAK HAOO 
FKHOYOMAPMOM TAM HHMNHANDEHSTKHNONDN FTF OCOD FRAC HDANPOPM ANDHOOM OWN 
DEE MNMMNAOUIE ANNO & FOC ANE HANDOVO ST ADK OME NODAMNMNAIEN AN ANADO 
OFAMNOMNOPM AMO PHAPNAINAN&S SHOMAFDNANAOMNOHOMMOarAPONNSOSFOOMA 
BDNMNNYNOODO sreerNtSFCOMDIODSOOQCODOLFOMHENVOODADAPOOCTS NOM Ort 
DFDOP MM MDD IONNA POM FOOD TOF NDOE OW) FMNMYAVET Me SOCK NN Sst 
eeoeeeveevreevrees eee teespeveee* ert eeevreseeveeeeseereeveeeteeseeentee 
CTRL RL ne ee ee ee eT 
' ( ! os i} t ' 
DOM MME NNN NG FOTO HO MRNA ANNAN NAAANAMMMMMMMAMMNMaAMMmAnM 
oleolelelelelelelelelelelelelelelel(elel ela ele eloleleleleleleoeleolelolelelslelslelelolelelelele) 
+H Arr rere errr errerereereerrrrtrerrtrrer)erteerbeete tere) ee 
ARARAARARAARAQANgQAAARANAAARARAAAAARAAAAARAAAAAAARAAaANAANRAARARAARDa 
CANDAMOMN GS FONE ATK DOM SFOONON ANE SOMMNANNONS MODOO Fo TOR OOO 
DNNNANDT ONNIDAMN OMS AMOM DOK POMONNOHNOIDHNOAOREKAUMOCOMemON 
DPONANNADS FAOM™ KODMEH ODOM KH HOOWMIM AD SAN INEMOKEMNOODNODDO OM 
DOANNDOCOM OFOME OF OD FOMNANANANT NT DOM ANMOMM FHA UM HODHANIM 
OPO FMM (NOQOMNDON F DINWM OF OM FE MOMNOMDINIIDATNOMEM IME OOUFODNATH WH 
CDAD AQOWUINATNDO YOM IF MNO DUIMNT &NNDBMMIMOUMN OTK NODA 
PDO DSTA OS FDS TFG H KOO KKK NOM OM AFMNNAT HOPE ETMMNOOF EE OWUMWM 
eeoeoee#eerteertestpeep ®eeerterset eee tees teeetrerteteteeeestste ft eeeeeweee 
So eo TTT TEP eT MITTEN (rte OT 
' 


leleleleleleleleleleleleleleleleleolelelelelelelelelelalelalelelelelelealelelele slolelelelelelala) 
e*esv8esstpeeoe? fF eeeerteeeestereeeseeeereerteeeee#e*eeeeefeef# fe oes * @ ef @ @ & @ 


lelelelelelealeleleleleleleleololelelelel ele) cle elelelelelelele elelelelslelelelelelselelelelelele) 


MNAFNOE DOTNMAOM VHONM FWOWDHNOSCNAWOPVOTAMAOP VDHNONM=ANO 
DO OO OOOO FSSC ECP KTANNANANANANMMM MM OM 0) =F xt at SPT a ot SCN 
Lr er 2 ee | 


ielelelelelelelelelelelelelelelelelelslelalelelelolelolelelololelelcleloleilcleloleleleleleleololo) 


94 


Ge See eS 
. 


Oooo 


5 
ae) 
-6 


owm 


Soke go Gk took kk ok ok to dk tok ok ac RR dk hk 
If YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 


-TYPE *Y¥* TO EXIT- 


(ANY OTHER INPUT TO CONTINUE) 
JR look dak ok kia ok kkk fk tok gtok tok kook 


B (120) 5/0 
C {131 R/O 
E (122) 5/0 


--- Your Fortran program _is now being loaded ... 
ees Execution will soon follow ... 
EXECUTICN BEGINS... 


THIS 


oO 
FO Uy GI 0 


L INPUTS, 
L Ts, 
ULATED. 


TO REVIEW 
AVED AS DAT 


wo 


are 
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CLEAR THE SCREEN TO CONTINUE. 
THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. PLOT THE DATA YOU JUST CALCULAT 
2. PLCT A CURVE THAT YOU PREVIOUSL 


ENTER 1 OR 2 


ED. 
Y SAVED. 


5 
; 
YCU MAY PLOT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DO YOU WISH TO PLOT? 
3 
WHICH TYEE CF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMEER 1? 
1. STATE VARIABLE (IE., X1, X2, ETC) 
2: FEEDBACK CCNTRCL (If., 6 = -c¥*x) 
3. CONTROL INPUT (IE., Ul, U2, ETC. 
4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAT2, ETC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC) 
ENTER 1,2,3,4 OR 5 
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AD-A15S2 224 DEVELOPMENT OF GRAPHICAL TIME RESPONSE “USING THE 
op eta Lt ert NAYAL POSTGRADUATE SCHOOL MONTEREY 
UNCLASSIFIED F/G 972 


egies 
—_ 3s 4 
= 
1.2 

2S WA Mes 


MICROCOPY RESOLUTION TEST CHART 
NATIONAL BUREAU OF STANDARDS 1062 A 


WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YCU WANT TC PLOT AS THE NUMBER 1 CURVE VS TIME? 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NOTE: 1. 4O ee eecaaee oe LENGTH 
2. GREEK SY LS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED i PARENTHESES. 
IE. (A) => ALPHA 
E) => BET 
F) => PHI 
C) => THETA 
angular position - {(x)1 


WHICH TYEE CF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 
J. STATE VARIABLE rise, X2y at 
2: REEDBACK CONTRG 
- CONTROL | (IE Ce) 
> STATE MATE iatle {IELs HATA KHAT 2, ETC.) 
. STATE ONSTRUCTION ERROR (1E., X1-XfaT1 
Sihts EEC ETC) 
: ENTER 1,2,3,4 OR 5 
2 
WHAT IS THE SUBSCRIPT OF THE STATE ESTIMATE THAT 
YCU WANT TC PLOT AS THE NUMBER 2 CURVE VS TIME? 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NCTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
B) => BETA 
F) => PHI 
2 Q) => THETA 
angular position estimate - (x)e1 
WHICH TYEE CF VARIABLE DO YOU WISH TO PLOT AS CUSVE NUMBER 3? 
1. STATE VARIABLE (IE., X1 re ETC) 
2. FEEDBACK CONTRCIL (zE. 
3. CCNTRCL INPUT (IE. ot. 2, ET Ce) 
4. STATE ESTIMATE dopser ve 8) (TEs HAT XHAI 2, ETC.) 
5: STATE RECONSTRUCTION ERROR , Xi-xfhatt, 
X2-XHAT2, ETC) 
ENTER 1,2,2,4 OR 5 
2 


WHAT IS THE SUBSCRIFT OF THE FEFDBACK CONTROL THAT 
YOU WANT TO PLOT AS THE NUMBER 3 CURVE VS TIME? 
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WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NOTE: 1. 40 CHARACTEES Le LENGTH 
2. GREEK SYMBOLS and BE Eeryree FOR ANY Lz2TTEKS 
ENCLOSED IN DARENT ESES 
ITE. (A) => ALPHA 
B) => BETA 
F) => PHI 
Q) => THETA 


input voltage - v 


YCU MAY USE UP TO 3 HEADINGS. 
HCW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 


0, 1, 2 OR 3 


WHAT IS THE DESIRED HEADING NUMBER 1? 


NCTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE soe FOR ANY LETTEES 
ENCLOSELC IN PARENTHES 
IE. (A) => ALPHA 
B) => BETA 
F) => PHI 
Q) =>_THETA 


filter + regulator closed loop system 
WHAT IS THE DESIRED HEADING NUMBER 2? 


NCTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WIIL Be PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALFHA 
B) => BETA 
F) => PHI 
Q) => THETA 


example 5.1 
WHAT IS THE DESIRED HEADING NUMBER 3? 


NCTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
B) => BETA 
F) => PHI 
Q) => THETA 
linear oftinal contrc¢el syst 
>> USING A PRE-ALLOCATED DATASET FORK UNIT FT17F001. 
THE FOLLOWING OFTIONS ARE AVAILABLE. 
1. BEGIN NEW GRAPH OF CTHER CCNTROLS, STATES, OR ESTIMATES. 
2. REFLOT PREVIOUSLY SAVED GRAFH DATS. 
3. ELIT THE CURRENT GRAPH 
4. PLOT REVISED GRAPH ON THE TEK618. ; 
5. gort ANC /OR era METAFILE OF THE CURVES. 
REVICUSLY SAVED 
ELECT A NUMBER BETWEEN 1 AND 5S. 


THE GRAPH EDIT MENU 
1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT FLO‘. ee 
22 DELETE CURVE FROM CURRENT PLOT. py 1g 
3. ELIT CURVE TITLE(S “ine 
4. EDIT PAGE HEALING (S) « one 
5. CEANGE THE Y-AXIS SCALE. : 
6. CHANGE THE TIME AXIS SCALE. a 
7. CHANGE PLOT SIZE. (DEFAOLT 1 IS 8.5 X 6.9) ; 
82 CHANCE THE LETTERI IGH 
9: CHANGE POSITICN OF THE LEGEND. 
10. EDITING COMPLETE. s 
SELECT A NUMBER EETWEEN 1 AND 10. 
9 
HOW MANY INCHES IN on Be y%y PERESTION 
LEFT Of RIGHT), D YOU WANT TO MOVE 
OVE THE LEGEND BCX FROM ITS PRESENT POSITION . < 
NOTE: 1. DEFAULT PLOT SIZE IS &.5 X 6.0 
2; LEFT IS NEGATIVE 
: 3. RIGHT IS POSITIVE 
0 
HOW MANY INCHES IN THE Y DIRECTION «a 
fue OR DOW) DO YOU WANT TO MOVE Ph 
OVE THE LEGEND BOX FROM ITS PRESENT POSITION 
NOTE: 1. DEFAULT PAGE SIZE IS &.5 X 6.0 
> DOBN IS NEGATIVE 
: 3. UE IS POSITIVE 
y; —- 
THE GRAPH EDIT MENU ae? 
1. CHANGE VARIABLES OR ADD A CURVE ON THE CUREENT PLOT. 
22 DELETE CURVE FROM CURRENT PLOT. 
3: ELIT CURVE TI ILE (S) 
4. EDIT PAGE HEADING ( he a 
5: CHANGE THE Y-AXIS SCALE. ele 
6. CHANGE THE TIME AXIS SCALE. 
7. CHANGE PLOT SIZE. (QEEAOLT IS 8.5 X 6.0), 
8. CHANGE THE LETIERING HEIGHT. 
9. CHANGE POSITICN OF THE LEGEND. 
10. EDITING COMPLETE. 
SELECT A NUMBER BETWEEN 1 AND 10. 
5 
ON WHICH CURVE CO YOU WANT TO CHANGE THE Y- SCALE? 
. ENTER CURVE NUMBER- 1, 2, 3, OR 4 
; 
- WHAT IS THE NEW Y-MIN VALUE AT THE ORIGIN? 
-.025 
7 WHAT IS THE NEW Y-MAX VALUE? 
0.125 


THE GKAPH EDIT MENU ee 


$8 


OR ADD A CURVE ON THE CURRENT FICT. 
CO RRENT PLOT. 


HQ 
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PQ 
MWHWHHAWO <3 
NOmte wwe dope 


XIS SCALE. 
: {REA OLT IS 8.5 X 6.0), 
ie G HEIGHT 


3: F THE LEGEND. 
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SELECT NUMBER BETWEEN 1 AND 10. 


10 
E FOLLOWING CPTIONS ARE AVAILABLE. 


EW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
REV AVED GRAPH DATS. 


1. 
2. 
3. APH. 

4, ON THE TEK618. 
5. ETAFILE OF THE CURVES. 


OmreHt 


SELECT A NUMBER BETWEEN 1 AND 5. 
4 


a 
= 


<SHRHIH 


G OFTIONS ARE AVAILABLE. 


ER CONTROLS, STATES, OR ESTIMATES. 
D GRAPH DaTv&. 

E TEK618. 

E OF THE CURVES. 


= 
e 
PRO bi ht bod 
watt ttim 
BHOHMA 
2 UObrrhit 


SELECT A NUMBER BETWEEN 1 AND 5. 


5 


WHAT FILE NAME [C YOU WANT THE CURVE DATA STORED UNDER? 
: (8 CHARACTERS MAX) 
filtereg 


THE CURVE [ATA IS BEING FILED UNDER FILTEREG DATA 
END OF DISSELA 9.0 -- 26332 VECTORS GENERATED IN 2 PLOT FRAMES 


FROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 
5020 VIRTUAL STORAGE REFERENCES; 9 SEADS; 0 WRITES. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 
2. QUIT. 
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' 
nth Mme 


to neat 


ENTER 1 OR 2 


WHAT FILE NAME IS THE DATA STORED UNDER? 


filtereg 
THE CURVE DATA IS BEING LOADED FROM FILE FILTERES DATA 2 
>> USING A PRE-ALLOCATED DATASET FOR UNIT FT18F0)1. fe 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 
2. QUIT. 


ENTER 7 OR 2 


32 GENE 
E F ISSCO, S 
FE 9 READS; 0 


ro 


ORS EATE 
Q ’ 
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HHH HH HK RK KK KK KK RK KKH KKK EK KKK RK KK RK KKK KKK KKK KK KK Ke 


DO _ YOU WANT A VRSTEC PLOITER SMOOTH COPY OF THE 
TEE DISSPLA eee ATT oR Ny YOU JUST CREATED? 


He HK He ee He ee He HH Me he eke eK a KK RK ee a ok ea oe ee tok tok & 


120) 4/9 “ties 
DASD 001 LINKED R/O; R/W BY MVS; R/O BY 0700P era 
Z §99 4} R/C - OS Tle 
DASD 001 DETACHED tye 
CREATING NEW FILE: 
CREATING NEW FILE: c 
PUN FILE 8317 TO MVS COPY 001 NOHOLD 

DASD 120 DETACHED 


PFE PEE SESE SELES SS TR SESS ELI SS TSE E SST ST TS ESE PS STE TEL RS ESE SELES ESS 


YOUR GRAPH(S) CAN BE PICKED UP AT THE COMPUTER CENTER. 


TEE GRAPH (S) WILL BE ADDRESSED TO "POP (USER ID)". 
HHA RARE EEE BK KK RK AK EERE ERK E EE EEK EERE RK RK RK EE RR 


ERE RARER ARK EK RE EKER EK RK ERK REE EK KK RK RK KE KKH HR KK KK 


DO YOU WANT TO oe 3 


1. RUN OPTSYSX AGAIN 
2. RUN THE ELOT PROGRAM US 
fer ANOTHER TYPE O 


ENTER 1, 2 OR 3 
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3 
‘6 PES ERE ESTEE SEES SSS TES ETT EP ES ST TESTS TS TESS ESS STS STR SST ESE LETS SS : 
HAVE A GOOD DAY!! - 
= PEPE SSE LESS ETT S ES REET ES ECS STE E SES ESP SESE TSS TST STS SEES EET SS SS IE 
| R; T=19.00/31.53 09:37:38 


record of 
END RECCKRDING OF TERMINAL SESSION, 


The graphical output generated by this exanple follows 
as figure 3.4. 
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A. CCNCIUSIONS 


As an ultimate evaluation of the computational abilities 
of OFISYSX, the program was tested usiny an 82 X 82 matrix 
of aircraft longitudinal motion equations for the %-29A 
experirental forward-swept wing Fighter aircraft prototyre, 
provided by NASA-Edwards. 

For a system of equations of this magnitude, all program 
arrays were re-dimensioned, and a 2-megabyte virtual machine 
size was required. 

The graphical time response curves generated from the 
X-29A longitudinal system matrix follow as figures 4.1, 4.2, 
4.3 and 4.4. The accuracy of these time response curves is 
mixed. All of the states shown have the correct waveforms, 
but differ ina scale factor of approximately times 10.0. 
Unfortunately the data supplied by NASA was not explicit 
regarding hew the control input was applied, and whether any 
additional gains were used in their simulation of the 
system. Time constrains did not allow the clarification of 
these prcblem areas, however the results of the X-29A longi- 
tudinal system are encouraging and should be a topic of 
further research. 
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It is hoped that ccntrol system instructors will 
encouraye their students to use this interactive graphical 
time response pregram for all applicable class projects; and 
that its enhanced capabilities will stimulate both interest 
in and research on basic systems control problems, as well 
as more advanced designs. 


EB. RECCEMENDATIONS 


Based on the results of this thesis, three areas emerged 


as possibilities for further research and study: 
1. Exrogram Availability 


The use of OPTSYSX and similar design programs 
should be encouraged in all undergraduate and graduate level 
courses involved in the analysis and design of centrol 
systems. Toward this end, it is recommended that the OFTSYS 
time response programs be placed in the non-IMSL likrary of 
subroutines, making it easily available to all pctertial 


users. 
2. FErogram Memory Requitements 


When configured for large matrix operations (98X98) 
the CFTSYSX program requires over 2 megabytes of virtual 
memory. Virtual machines of this size are not normally 
available te a user. The memory usage for matrix storage is 
a possikle area of improvement in the efficiency of the 
CPTSYSX program design. All matrices calculated in  TSYSX 
{except DUMMY matrices) are still available when the run is 
finished. This Simplifies program operation but uses an 
excessive amount of memory. Memory usage should re studied 
and frogram modifications should be made to reduce the 
excessive memory requirenents. 
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3. Further 


Modifications 


Program sequencing during optimal filt2r synthesis 
should be modified. At the present time a regulator must be 
designed or supplied when a filter is designed. Various 
test runs indicate that this problem can be overridden if 
the number cf controls (Nc) is given as zero, but this is 
not a viakle soluticn for systems which use a driving 
function. 


109 


APPENDIX A. OPTSYS EXEC LISTING 
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THE OFTSYS EXEC 


CONTROLS THE OPTSYSX, OPTCALC AND OPTIPLOT 
TC DETERMINE THE TIMf RESPONSE OF A SYSTE 


BY H. A. DIEL 
VERSION 1.0 16 JULY 1934 


eR eo eK KK RRR KKK RK KKK KKK KK KKK KEK REE KK RK EK KKK 
Pit ESE SESE EET CSS ETE SES ES STS ETS ST ESS SPE SSP PES TE SESS TE EEE LSS 


RRS HEHEHE 


* * 
* CHECK FOR USER'S VM SIZE = > THAN 1 MEGBYTE : 
HR RRR I FO RI IER a i tok doe 
VMSIZE 

SIF &FC GE 1024 §&GOTC -TWO 


CLESCEN 
EBEGIYPE -~ENDTHREE 


WH HE He eH He He He eH He oe A he He Re he ee he ea oo oe eo eo ee oe a eo ee tke & 


YOU MUST HAVE A 1M CR LARGER VIRTUAL MACHINE 
TO RUN THIS OPTSYS PROGRAM 


TO DEFINE A 1M VIRTUAL MACHINE: 
DEFINE STORAGE 1m 
I CMS 
OPTSYS 
FOR SYSTEMS LARGER THA 
OETAIN A LISTING OF THE 0 
AND FOLLOW INSTRUC 


HERA ADEA EAR RK EK RRR ERE RE KEK EK 
~ENDTHREE 
SEXIT &RC 

WO 


CLRSCRN 
EBEGTYPE -ENDZERC 


THE CPISYS EXEC CCNTROLS A TRIO OF PROGRAMS: 

1. OFTSYSX FORTEAN SYSTEM ANALYSIS 

2- OPTCALC FORTRAN CALCULATE TIME Berra 

3. OFTPLOT FORTEAN DISSFLA PLOTTING ROUTINE) 
EACH FRCGRAM PASSES INFORMATION TO THE NEXT PROGRAM _ 
THROUGH A DATA FILE WRITIEN TO THE USERS DISK. IN THIS 
CASE, THESE FILES ARE "OETMAT DATA" AND “OPTPLOT DATA". 
THE SIZE OF THESE FILES VARY WITH THE SYSTEM ORDER, ANE 
CAN USE ABOUT 20% CF THE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 


Ree ee ea ies Eg re me Sar ae ae ge Site SQe SS meer tee og ay 
fa Oe et Ere Te ee. Bese e Mas Sig TS a) Sette eee eral eee tee 


Pee ech Sey eae Sts Wh ee Tei ae Te 


rd 


- TYPE "E" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 

-ENDZERO 

SREAD VAES_§ANS - 
SIF .§ANS EC .E SEXIT ERC a 
RRRK KE RKKGKK KKK KE KKKKRKEBREKEKK EKER KEKE KKK EKEKKEKEEK ERK EES r) 
* * 

* ALLOW THE USE OF AN OLD "OFTMAT DATA A1" * 


* * 

RAR EES HR RE EE ARK EERE AE EERE EERE RR EERE EE RE KK 

RENAME OFTMAT pei Al peaore DATA Al 

SIF &EC KE 0 &GO ~FIRS 

pee OFTSYS BATA Al OpTMAT DATA Al ; 
EBEGIYPE -ENDONE 


PES SESE SSIES SETS SET ESE SESE TES ESS TESTE SSE TET ESTES TTT RESTS SEES TS 


You stHad DATA FILE NAMED ‘OPTMAT DATA* ON YOUR 


DISK T WAS pee aoorae GENERATED BY THE preys 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT FUN. 

IF_ You youRe LIKE T0 USE THESE SAME MATRICES FOR 
THIS KUN, THE OPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROPRIATE TIME, 


IF You TYPE (Y) ES. ’ 


ANY OTEER INPUT WILL RESULT IN THAT FILE BEING ERASED! 
SHH HR RRR RR FR RR RK RRR RIO RIOR BR FR TR I 


-EN DONE . 
SREAD VARS ANS 

BIF GANS EC .¥ &GOTC -ONE > 
SEI OR RR at ik tk i RoR tk Rk Rk HR RK RK i tk EK : 

* * 

* ERASE THE OLD "OPIMAT DATA AI" DATA FILE * 

*# PLACE "000 0" IN TEE NEW "OPTMAT DATA FILE" * 

* TO ACT AS A FLAG FCR OPTSYSX AND OPTCAILC * 

PEE ESSERE ESS ST ERE ESE SESE SS PT RSS EES ESLER SEES SS ELST LESS LS TEES SYS ) va 
ERASE OPTMAT DATA A1 

ESTACK 000 0 

FILESICK OPTMAT DATA A1 P 80 1 

-THIRD 

CLRSCEN 

SBEGTYPE -ENDFOUR 


WAH HK KK HR KK KKK ERE KK KER KKK KKK ERE KEKE KEKE KER EEK EEK 


DC YCU WANT THE NUMERICAL OUTPUT FROM CPTSYSX TO GO i] 
TC YOURS TERMINAL 1o GR DS CR TO A D(ISK) FILE? : 


RHEAKKEAKKHEAK KKK KK KE KAKKEKKK EEK KEKE EERE KKK ERE KEK KEKEKEKK ESE ea ae 
-END oe 


EEEAD VARS €ANS 
SIF .§ANS EQ .S &GOTC -FOURTH . ’ 
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we a oe ce Oe ae eget an ie 


NS EC .D &GOTC -FIFTH 
E ~ENDFIVE 


PES ESEET EEE ETO RES ESSER SS TESTE SESE TES ETE SES STS CESS EEE TEES SS FS 


YOU MUST ANSWER S(CREEN) OR D (ISK). 
Prt SS EES ERLE STS SETS SERS ESTES ESE SS ES SESS ST TCS ESS TST SS ESET SE EES tt 


~ENDFIVE 
CP SLEEP 3 eae 


EBEGTYPE -ENDSIX 


PILE SESE LES ETO SE EOS SETI ETI T STP eT errr re Per Te TT eT Te eet ee eres SS 
OUTPUT WILL COME TO YOUR TERMINAL SCREEN. 


BRK RR RO tO dk tok fk I OR dO a do 2k 

~ENDSIX 

CP SLEEP 1 SEC Y 
ETYPE LOADING OPISYS.... . 
FILEDEF 06 TERM seeree FA BLKSIZE_ 133 tee 
FILEDEF 8 DISK OPTPLOT DATA A1 genes ee 
FILEDEF S DISK OPTMAT DATA Al (PERM re 
GLOBAL TXILIB FORTMOD2 MOD2EEH IMSLDP NONIMSL oe 
LOAD OPTSYSX (START uke 
&GOTC -FIVE ; 
-PIFI# 

CLRSCEN 

EBEGTYPE -ENDSEVEN 


RHR KEKE EKHRK RKK KK KK KE RAE ERK ERK KEKHE RE KEKE HERE EKKK ERKEK SEEK EK 


CUTPUT WILL GO IC DISK FILE ‘OUTPUTX LISTING Ai'* 
BRI KR RK BRK KI RRR RR FRR RO FOR IR FOR BOR SOR OK tO 


~ENDSEVEN 

CP SLEEP 1 SEC 

ETYPE LCADING OPTSYS.... 

FILEDEF C6 DISK OUTPUIX LISTING A1 

FILEDEF 8 DISK OPTPLCT DATA Ai Agaes 

FILEDEF 9 DISK OPTMAT DATA Al (PERM 

GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP NONIMSL 

cea acre (START = 


E ~ENDEIGHT 


FESS SESE EEC ELIS SESE SESE ES SSSR SSE SSE TEP STE TS SE TSS See SES TES ESS 


IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 


“TYPE "Y¥*' TO EXIT- 


ANY OTHER INPUT TO CONTI NUE) 
PET ESSE SSO SII Se SP er eS Te eT Pee STP Te eT TET Tete ee tee ee SS 


-ENDEIG 
GSREAD VARS 

eIF gGANs Ee NY SEXIT ERC 

RAK RH KK RK KK KKK KEE EK EK KKK KKK KK KK KK ERK KEKE KKK KE RT HK 
* * 
* CRECK FOR DATA IN THE FILE "OPTMAT DATA " * 
* BEFCRE LOADING OPTCALC * 
WKH KH RK EEK KEK KK KK KK ARK HRKKEK RHEE RK KEKE KEKE KKK REE KKH EK EK HK 
FINDSTAK OPIMAT DATA A1 CO1 0 LIMOO2 ALL GROUP1 

SREAD VARS &A1 &A2 

SIF .&A2 EQ .0 SEXIT ERC 

STYPE LCADING OPTCALC.... 

FILEDER § DISK OPTPLCT DATA a1 (PERM 

FILEDFF 9 DISK OPTMAT DATA A1 (PERM 

GLOBAL TXTLIB FORTMODZ MOD2EEH IMSLDP NONIMSL 

LOAD CETCALC (START 

~SEVENTH 

CLRSCEN 

SBEGTYPE -ENDNINE 


TEP SRE SE SESE LES SET EEE SESE TE SESE TEE ES TELE STS ETS TELS SES EE SS tS 


IF YOU ARE DISSATISFIED WITH THE RESULTS 
THOS FAR AND WOULD LIKE TO EXIT TO CMS, 


“TYPE *Y* To EXIT- 


(ANY OTHER INPUT TO CONTINUE) 


RAAB AR REE AA AK RAE RK RRR ER EE RK EE EO KE EE 
~ENDNINE 

EREAD VAES &ANS 

cIP_ .é&A 


e SANS eee oY SEXIT &RC 
FILELDEF SK_OPTPLCT DATA A1 (PERM 
ESC BISSerA OPTFLOT 
HH EE ARK RAK ERE KEK EK RRR RR EER RRR ERE RK FE ROE OK EEK 
* * 
* CHECK FOR FILE "DISSPLA METAFILE A&" ON * 
: THE USER'S CISK BEFORE GOING TO DISSPOP . 


OO PETET ESET ET EST TTTETOTTETETETEST TTT C TT TEST CTT ee eT ere TT Te 
RENAME Dee rLe METAFILE - OPTSYS METAFILE A4 

SIF &RC 0 &GOTO -TENT 

RENAME oftsY¥s METAFILE Ab DISSPLA METAFILE A4¥ 
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ee ete 


E -~ENDTEN 


TET eCCCETOCT SOC SeS SECT OC STL SLT CCC CLS Te STS Tee Te TTS See SE ee TS eee 


DC YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE : 
THE DISSPLA Tey OR a YOU JUST CREATED? ' 


REA KK EAR AEE KE RK KK KE GEE KKK KK REREEKRKEK KE RK EN KKK KEKE ER KKK HK EK 


@ ' 
PPS SeL EES ES ELT ESS LSS EELS ESE RSE SE SES SL SESE SESE TESTIS SESE SEES ES 

: YOU MUST ANSWER Y(ES) OR N(O0). 
RIK KKH RH RR KK KK KKK KKK KK KKK KHER KK KEK KK KE KEKE ERK ERE fos 
-ENDELEVEN Bin 
CP SLEEP 4 SEC ees 
&GOTC -EIGHTH ae 
-NINTH Oe 
EXEC DISSPOF VRSTEC cae 
LRSCRN ae 
EBEGTYEE -ENDTWELVE a 
FESS SESE EEE ETT ES ECE EE ESS TTS SESS SST SESS ESTES ST SLT STS SES TEES SSS : 

YOUR GRAPH(S) CAN EE PICKED UP AT THE COMPUTER CENTER. 
THE GRAPH(S) WIIL BE ADDRESSED TO "POP (USER ID)". 

RK KKH KKK KKH KKK KK KK KKK KK RK KEK KKK KE KKK KK RAKE KEKEK EEK EKK 

-END TWELVE 

le CP SLEEP 5 SEC . 
=TENTH 
CLRSCEN 
SBEGTYPS ~ENDTHIRTEEN 


PES SEE SEES ES ETE EES ESE SESE TEESE ESE LES TESS SLES TET EER S SSE SS TS SS 
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DO YOU WANT TO 


1. RUN DPTSYSX i 
2: RUN THE PLOT PROGRAM USING THE SAME MATRICES? 
C PLOT ANOTHER TYPE OF SYSTEM (OPEN/CLOSED) ) - 
3. QUIT 
oa ENTER 1, 2 OR 3 
ee PESTS SSE SELES ES SSS TESS ST SS SESS SESE ESE SEES TS SS UTS TS STR SES SS SS TS ee 
me -END THIRTEEN pas 
7 SREAD VAES &ANS 
BIE GANS EC .] §GO70 -THIRD 
SIF [ANS E€ 22 §GOTC -SIXTH 
CLRSCRN 
; SEEGTYPE -ENDGOODBY 


WRK Fe HH HK HR KR RR KK RK RK RK KKK HK EK KE 


HAVE A GOOD DAY!! 


We HCH He he He He HH he ee eK eH KK KK KK KK RK KK KK ke & 


~ENDGCODEY 
cP Sr Eee 3 SEC 


N 
SEXIT ERC 


ow 
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